Debye Sources and the Numerical Solution 
of the Time Harmonic Maxwell Equations 



Charles L. Epstein* and Leslie Greengard^ 



January 14, 2008 



Abstract 



In this paper, we develop a new representation for outgoing solutions to 
the time harmonic Maxwell equations in unbounded domains in M 3 . This rep- 
resentation leads to a Fredholm integral equation of the second kind for solv- 
ing the problem of scattering from a perfect conductor, which does not suffer 
from spurious resonances or low frequency breakdown, although it requires 
the inversion of the scalar surface Laplacian on the domain boundary. In the 
course of our analysis, we give a new proof of the existence of non-trivial 
families of time harmonic solutions with vanishing normal components that 
arise when the boundary of the domain is not simply connected. We refer to 
these as /c-Neumann fields, since they generalize, to non-zero wave numbers, 
the classical harmonic Neumann fields. The existence of /c-harmonic fields 
was established earlier by Kress. 

Introduction 

Electromagnetic wave propagation in a uniform, nonconducting, isotropic medium 
in M 3 is described by the Maxwell equations 
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V -£(x,i) = 0, V-H(x,t) = 0, 

where £ and TL denote the electric and magnetic fields, respectively and e, fi are the 
electrical permittivity and magnetic permeability of the medium. We will restrict 
our attention to the time-harmonic case and write 

£(x,t)=Rel^^-e-^ t ) and H(x,t) = Re( ^M e -^j , (2) 



V~e J I ^ 

The superscript is used to emphasize that E tot and H tot define the total electric 
and magnetic fields, respectively. In electromagnetic scattering, they are generally 
written as a sum 

E tot (x) = E in (x) + E(x), H tot (x) = H in (x) + H(x), (3) 

where {E m , H m } describe a known incident field and {E, H} denote the scat- 
tered field of interest. With the scaling in ©, the Maxwell equations take the 
simpler form 

V x H tot = -ikE tot (4) 

V x E tot = ikH tot . 

We are particularly interested in the problem of scattering from a perfect conductor 
in an exterior region, which we denote by Q,. For a perfect conductor |[T3ll22l . the 
conditions to be enforced on T, the boundary of Q, are 

n x E tot = 0, n ■ H tot = 0. (5) 

The scattered field is assumed to satisfy the Silver-Miiller radiation condition: 

lim (Hx --e) =o( - ) . (6) 

r^oo \ r ) \r J 

This problem has been studied rather intensively for many decades, and we do 
not seek to review the literature here, except to observe that there are two distinct 
approaches in widespread use. When the scatterer is a sphere, a simple and elegant 
theory exists due to Lorenz, Debye and Mie (6l[T0l[T6l[T9l- It is based on two scalar 
potentials (generally called Debye potentials), and the mathematical machinery of 
vector spherical harmonics. In particular, one represents E, H as 



E = V x V x (rv) + ikV x (ru) 

H = V x V x (ru) - ikV x (rv) (7) 
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where the Debye potentials u, v satisfy the scalar Helmholtz equation 



An + k 2 u = 0, Av + k 2 v = , 

with Helmholtz parameter (wave number) k 2 = u> 2 e/x. Salient features of this 
approach are (a) that the boundary value problem 

n x E = -nx E in (8) 
n H = -n H in (9) 

is uniquely solvable for any k with non-negative imaginary part and (b) that, as 
k — > (uj — » 0), the electric and magnetic fields uncouple gracefully. In the static 
limit, E is due to the scalar potential v alone, which is, in turn, determined by the 
boundary data —nx E m . Likewise, H is due to the scalar potential u alone, which 
is determined by the boundary data — n ■ H m . 

For regions of arbitrary shape, on the other hand, integral formulations of the 
Maxwell equations are generally based on the classical vector and scalar potentials 
(in the Lorenz gauge): 

E = ikA-V(/> (10) 
H = V x A (11) 



where 

<j>{x 

with 



A(x) = J g k (x-y)J(y)dA y 
) = ^ J^g k (x - y)(Vr ■ J)(y)dA l 



e ik\x\ 

9k(x) = —— 
4tt x 



Here, J is a surface current (a tangential vector field) and Vr • J denotes its surface 
divergence. 

Maue [17'] proposed the electric field integral equation (EFIE) for the unknown 
current J by enforcing the condition ([8]> using the representation ([TOl . Because of 
the V4> term, however, the result is a hypersingular equation. Maue also proposed 
the magnetic field integral equation (MFIE), based on ((TTb . The boundary condi- 
tion for H can be derived from the Maxwell equations and an appropriate limiting 
process on the surface of a perfect conductor lfT3l l22l : 

J = nxH in + nxH, (12) 
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where n points into 0. Enforcing this condition for the unknown current J yields 
the MFIE, a second kind Fredholm equation. Unfortunately, both the MFIE and the 
EFIE have spurious resonances; that is, there exists a countable set of frequencies 
{kj} C R for which the integral equations are not invertible. As the {kj} are 
the eigenvalues of a self adjoint, elliptic boundary value problem on the bounded 
complement of Q, they are often referred to as interior resonances. Below the 
smallest such kj, the MFIE is well-conditioned. Spurious resonances, however, 
are only one difficulty. A second problem stems from the representation of the 
electric field itself. Unlike the Debye representation, the electric field does not 
uncouple naturally from the magnetic field as k — > 0. Note that in (fTOb . E involves 
one term of order k and one term of order k^ 1 . This results in what is referred to 
as "low-frequency breakdown" [36]. While low frequency breakdown is a more 
transparent problem in the context of the EFIE, the MFIE is not immune l35l . 
Knowing the current J is sufficient for computing H, but not the electric field. 
The normal component of E, for example, is determined by the electric charge: 

V r • J 

n-E = p=^—. (13) 

ik 

As k — > 0, accuracy degrades dramatically - a phenomenon called "catastrophic 
cancellation" in numerical analysis. 

This state of affairs is both odd and unsatisfactory. For the exterior of a sphere, 
there is a simple, clean representation of the solution based on two scalar unknowns 
that results in a diagonal linear system. It has no spurious resonances and, at zero 
frequency, decouples naturally (with no loss of precision) into magnetostatic and 
electrostatic problems. The standard integral equation approaches available for 
general geometries do not reduce to a Debye-like formalism when restricted to a 
sphere. Instead, a sequence of modifications have been introduced to address the 
three problems discussed above: the existence of spurious resonances, the lack of 
a second kind integral equation valid for all frequencies, and the loss of accuracy 
due to low-frequency breakdown. 

An important step in addressing the first problem was the introduction in the 
1970's of the combined field integral equation (CFIE) (20l|25l. The CFIE avoids 
spurious resonances by taking a complex linear combination of the EFIE and the 
MFIE, both of which involve the surface current as the unknown. It is not a Fred- 
holm equation of the second kind, however, and still suffers from low frequency 
breakdown. One alternative approach, due to Yaghjian [33] , involves augmenting 
the MFIE with the condition (© or the EFIE with the condition (fT3l . He showed 
that (for geometries other than the sphere) the augmented equations yield unique 
solutions at all frequencies. Of the many formulations that have been introduced 
to overcome spurious resonances, variants of the CFIE have emerged as the most 
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frequently used in practice. 

For the second problem, the principal issue is that of overcoming the hyper- 
singular behavior of the CFIE. One solution is to introduce electric charge p as an 
additional variable [27]. In this approach, one defines the scalar potential (j> by 



While the hypersingular term is avoided, one must solve a Fredholm integral equa- 
tion subject to a differential-algebraic constraint ([TBI) . During the last several years, 
several promising approaches have been developed based on the construction of 
preconditioners. Christiansen and Nedelec [7] designed effective strategies for the 
EFIE based on Calderon formulas and the Helmholtz decomposition. Adams and 
Contopanagos et al. (H made use of the fact that the EFIE operator serves as its 
own preconditioner; more precisely, the composition of the hypersingular operator 
with itself equals the sum of the identity operator and a compact operator. A com- 
bined field integral equation using this preconditioned EFIE is both resonance-free 
and takes the form of a Fredholm equation of the second kind. Preconditioners 
have also been designed through the use of high frequency asymptotics [2J. Unfor- 
tunately, the implementation of these schemes can be rather involved on arbitrary 
surfaces and, like the MFIE, they still suffers from a form of low-frequency break- 
down in the evaluation of E once the integral equation has been solved. 

Finally, the third problem - namely the low-frequency breakdown of the inte- 
gral equations - has generally been handled through the use of specialized basis 
functions in the discretization of the current, such as the "loop and tree" method 
of OT1I32I . This is a kind of discrete surface Helmholtz decomposition of J. As 
the frequency goes to zero, the irrotational and solenoidal discretization elements 
become uncoupled, avoiding the scaling problem that causes loss of precision. 

We have chosen to investigate a rather different line of thought, motivated 
largely by the desire to extend the Debye potentials to surfaces of arbitrary shape. 
In essence, the Lorenz-Debye-Mie approach is based on expanding the potentials 
u, v from © as 




(14) 



and imposes the continuity condition 



V • J = ikp. 



(15) 



n,m 
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where h n (x) is the spherical Hankel function of order n, and Y™{0,4>) is the 
usual spherical harmonic of order n and degree m. This separation of variables 
approach is clearly not suitable in general. From a mathematical viewpoint, it 
works because of the close connection between the Laplacian in M 3 and the surface 
Laplace-Beltrami operator on the sphere. It is also worth noting that the Lorenz- 
Debye-Mie approach is not equivalent to a Fredholm equation of the second kind. 
It is invertible, resonance free and behaves properly at low frequencies, but it is 
hypersingular. Numerical difficulties are avoided simply because it is in diagonal 
form. 

The features of the Debye potentials that we wish to retain are their symmetry 
and the fact that, at zero frequency, the system uncouples into separate electrostatic 
and magnetostatic problems. For symmetry, we begin by using both potentials 
(A,(f>) and "antipotentials" (A m ,<p m ) as a formal representation of the electro- 
magnetic fields l22l : 

E = ikA-Vcp-V x A m (16) 
H = V x A + ikA m - V^ m (17) 

where 



A(x) = 


J 9k(x 


- y)j(y)dA y , 


<t>(x) = 




- y)r(y)dA y , 


A m (x) = 


J^9k(x 


- y)m(y)dA y , 


(f) m (x) = 


J^9k(x 


- y)q(y)dA y , 



together with the continuity conditions 

Vr • j = ikr , Vr • m = ikq. (19) 

Such a symmetric formulation is commonly used for scattering from a dielec- 
tric. For the perfect conductor, it underlies the combined source integral equa- 
tion method (CSIE) |[T8l , where j and m are both assumed to be derived from a 
"parent" current distribution j: 

j = aj, m = (1 — a)n x j 

for some parameter a. More precisely, the CSIE is derived using ([TBI ) with the vec- 
tor unknown j and enforcing the condition ([8]). Like the CFIE, it is a resonance-free 
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but hypersingular equation. It is important to recognize that, in this construction, 
the unknowns j, r, m, q are no longer physical quantities, j and r correspond to 
fictitious electric current and electric charge, while m and q correspond to fictitious 
magnetic current and magnetic charge. Perfect conductors do not support the latter. 
If the "physical" current supported on the surface T is desired, it must be computed 
in a second step. From (fT2l) . for example, we have J = nx (H m + H). This will 
not be the unknown j. 

The second and critical feature of our method is that we will use r and q as 
unknowns and construct j and m from them in such a way that the continuity 
conditions ( fT9l) are automatically satisfied. In particular, for simply connected do- 
mains, we will let 

j = Vr^ + n x V r ^ m (20) 
m = n x j 

where 

A r ^ = Vp* = ikr (21) 
Ar* = V r ^ m = -ikq. 

We wil refer to Ar as the surface Laplacian or Laplace-Beltrami operator. (In 
geometry, this name is usually applied to — Ap so that it is a non-negative operator, 
but we will use the the convention above consistently.) In any case, we will obtain 
the Helmholtz decomposition of the currents on the surface by construction. This 
avoids the obvious cause of low-frequency breakdown, since we never compute 
the O(l) quantities r, q from the O(k) quantities j, m with its attendant loss of 
accuracy. 

An obvious drawback of our approach, of course, is that it will require the 
inversion of a partial differential equation on the surface of the scatterer to compute 
^ and ^ rn . It is interesting to note that Scharstein proposed an investigation along 
these lines some years ago [26], using only the electric current J, but a detailed 
investigation of the theory was not carried out. 

We show below that our representation yields a second-kind integral equation 
for r and q that, in the simply connected case, has a unique solution for all fre- 
quencies with non-negative imaginary part. Furthermore, it behaves gracefully in 
the low frequency limit, uncoupling into an electrostatic problem involving r and a 
magnetostatic problem involving q. Because of the connection with the Debye the- 
ory, we refer to r and q as generalized Debye sources. We also present an analysis 
of the (more delicate) multiply-connected case. 
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1 Geometric Preliminaries 



Definition 1. Let D denote a bounded (not necessarily connected) region in M 3 
and let £1 denote the unbounded component of M 3 \ D. We will refer to £1 as the 
exterior region and to its boundary as T. We assume, without loss of generality, 
that M 3 \ D has no bounded components (that is, holes within the interior of D). 

Using standard topological terminology, let us assume that D is multiply con- 
nected with genus g. Then there exist surfaces S\ , . . . , S g in D such that D \ 
U| =1 5j is simply connected and surfaces T%, . . . , T g in R 3 \ D such that M 3 \ D \ 
Uj =1 Tj is simply connected. We denote by Aj the boundary of Sj and by Bj the 
boundary of Tj (see Fig. [[])■ 

Remark 1. We will refer to the curves {Aj\j = 1, ... ,g} as A-cycles. (They 
form a basis for the first homology group of R 3 \ D.) We will refer to the curves 
{Bj \j = 1, . . . , g} as B-cycles. (They form a basis for the first homology group of 



we refer to it as having mean zero on that component. We denote by .Mr,o the set 
of pairs of functions (/, g) on T with mean zero on every component. 

Lemma 1 (The mean zero condition). Let (r, q) be generalized Debye sources 
defined on a boundary T. Then (r, q) £ Mr,o- 

This is proven in section l6\2l (see eqs. (l89l ) and (|90ll). 

In multiply connected domains, the Helmholtz decomposition (1201 ) is incom- 
plete. From Hodge theory, however, we can write a surface vector field j as an 
orthogonal decomposition (in L2), the Hodge-Helmholtz decomposition: 



j R = V r ^ + n x V r * m 

for some ^, where j H satisfies 

Vr • 3h = 0, V r • (n x j H ) = 0. 

See Appendix IA.3I and [29 ]. Such vector fields are called harmonic vector fields 
(dual to harmonic 1-forms). We let 



D.) 



Definition 2. Let Tj denote a component of the boundary T. If 




3 = 3r + 3h 



(22) 



m R = nxj R 



m H = n x j H 
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B-cycles 




A-cycles 



Figure 1: A multiply connected domain D, consisting of three components (two of 
genus 1 and one of genus 0), with exterior f2. Cutting along the surfaces bounded 
by the "A-cycle" from D makes it simply connected. Adding the surfaces Tj 
bounded by the B-cycles makes Q\D simply connected. 

and (as before) 



Harmonic vector fields arise, in essence, because the Laplace-Beltrami operator on 
vector fields 



has a non-trivial nullspace on multiply connected surfaces. The dimension of the 
nullspace of Ap is equal to twice the genus, g, of the surface. We may, therefore, 
choose harmonic vector fields {j l H \ l = 1, . . . , 2g}, which form an orthogonal 
basis,w.r.t. L 2 (T), for this nullspace. In the multiply connected case, we can then 
define the harmonic components of j and m by 



m = n x j . 



Ai 



■rJ = v rVr j-nx V r V r • (n x j), 



(23) 




(24) 



l=i 



m H = n x j H . 
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The space of harmonic vector fields will be denoted Tt 1 ^). 

Given that the Laplace-Beltrami operator is not invertible, one must be careful 
in defining ^ and \l/ m . From Hodge theory, however, we know that it is invertible 
as a map from the space of mean zero functions Mr,o to itself. We denote by Rq 
the partial inverse of Ap acting on this space. Thus, we replace (|2TT) with 

= ikR r (25) 
* m = -ikR q , 

where r, q are the generalized Debye sources. 
Example 1. Consider a torus in M 3 , parametrized by 

x(s,t) = [(R + r cos t) cos s,(R + r cos t) sin s, r sin t] , 

with the z-axis as the axis of symmetry. A straightforward calculation shows that 

•l _ 1 dx(s,t) 

3h ~ (R + rcost) 2 ds 

and 

3 2 H = n * 3h 

are both harmonic vector fields. Since the genus of a torus is 1 , they form a basis 
for the two-dimensional space of harmonic vector fields on the surface. 

Remark 2. Much of the formal analysis in this paper is simplified through the use 
of differential forms and homology theory. In order to be accessible to a broader 
audience, however, we state the main results using the notation of vector calculus 
and defer most proofs to Sections [6] and [71 where we do make use of the language 
and power of this theory. 

2 Uniqueness Theorems for Exterior Electromagnetic Fields 

Let us denote by C + the closed upper half plane: 

C+ ={z£C: Imz > 0}. 

Definition 3. A solution to the time harmonic Maxwell equations in 17 that satisfies 
the Silver-Muller radiation condition will be referred to as an outgoing solution. 

That an outgoing solution to THME(fc) is determined by either the tangential 
components of the electric or magnetic fields is classical (H : 
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Theorem 1. Suppose that (E, H) is an outgoing solution to the THME(k) in an 
exterior region £lfor nonzero k G C + . If either n x E or n x H vanishes on T, 
then the solution is identically zero in 0,. 

Although we will eventually address the problem of scattering from a perfect 
conductor (n x E = 0), we turn our attention for the moment to the Maxwell 
equations in exterior domains with normal components specified on the bound- 
ary. While this is not a standard physical boundary value problem, there is prior 
work on uniqueness and it is a natural starting point for the analysis of symmetric 
representations of the fields. 

Theorem 2. [Yee, 1970]. Let (E, H) be an outgoing solution to the THME(k) in 
an exterior region Qfor nonzero k G C + . Suppose T is simply connected, and that 



Then E = H = inVL. 

When the boundary has non-trivial topology, a rather subtle argument shows 
that, in general, this is not true. In particular, if the sum of the genera of the 
boundary components of the exterior domain is g > 0, then for all frequencies with 
non-negative imaginary part, there is a 2g-dimensional space of outgoing solutions 
to the THME with vanishing normal components. The existence of these fields was 
proven by Kress (Theorem [3] below). 

Remark 3. In the static (harmonic) case, this fact has been known for decades QUI . 
More precisely, at k = 0, the THME separate into the system 



solutions to which are called harmonic vector fields (if they decay at infinity). 
When their normal components vanish on the boundary, they are called harmonic 
Neumann fields. If their tangential components vanish, they are called harmonic 
Dirichlet fields. 

Theorem 3. [Kress, 1986 [151]. Suppose that (E, H) is an outgoing solution to 
the THME(k) in the exterior region Q,for nonzero k G C+. If every component of 
the boundary T is simply connected, then the solution is determined by the normal 
components n ■ E |r and n ■ H [r • If the sum of the genera of the components of 
r equals g > 0, then there is a subspace of outgoing solutions to THME(k) with 



n ■ E \ r = , n • H \ T = . 



V -E = 0, Vx£ = 



V-iT = 0, V xH = 



n ■ E |r= and n ■ H fr= 



(26) 



of dimension exactly 2g. 
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Lemma 2. [Kress, 1986 /T73T77. Let (E, H) be a solution to the THME(k) for 
nonzero k £ C + , in a region f2 with a smooth bounded boundary V. Then the 
normal components (n ■ E,n ■ H) lie in .Mr,o- 

Corollary 1. [Kress, 1986 [15]]. Suppose that (E, H) is an outgoing solution to 
the THME(k) with k £ C+ and satisfying 

n ■ E \y= f and n • H \r= h 

where f,h£ Mr,o are Holder continuous. Then (E, H) are uniquely determined 
subject to the specification of their circulations on the A-cycles (see Fig. [7]).- 

/ r • E ds = pj and / r • H ds = qj , (27) 

J Aj J Aj 

where pj ,qj € C are given numbers. 

Remark 4. We call solutions to THME(/c) that satisfy (l26l ) k-Neumann fields, and 
denote the space of such solutions by TLk{^l)- The conditions in Corollary Q] are 
familiar from the static (zero frequency) case, where g conditions must be specified 
for each of E and H separately, since the equations are uncoupled. For nonzero 
k, this symmetry is not required. We provide a different proof of existence in 
Theorem QTJ and a somewhat more general analysis of uniqueness in section loTTl 

Our representation also provides an effective means for numerically computing 
the /c-Neumann fields. These solutions will be needed in solving the problem of 
scattering from a mutiply-connected perfect conductor. 

First, however, we need to recall some classical facts about layer potentials. 



3 Jump Relations and Boundary Values of the Potentials 

In order to use the integral representations discussed above to solve boundary value 
problems, we need to find expressions for the restrictions of E and H to the bound- 
ary, in terms of the various potentials. In this section we collect the relevant results. 
Recall that T <^-> M 3 is a smooth, bounded surface (possibly disconnected). The 
unbounded component of M 3 \ T which we have denoted by $7, will be referred to 
as the "+" side of the boundary. The domain D (the union of the bounded com- 
ponents) will be referred to as the "— " side. We use n to denote the unit normal 
vector field along T, pointing into the unbounded component. 

The relevant limits are given in the following lemma, proofs of which can be 
found, for example, in (H. 
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Lemma 3. Let A and <p denote the vector and scalar potentials in ([78]) and let V 
be a smooth bounded surface in R 3 . For Xq 6 r, let x —> x^ indicate approach 
from Q. (+) or D (—), respectively, and let no denote the normal at Xq, with 
the normal derivative at Xq.. Then, for the scalar potential, we have 



1 



(28) 



where 



lirn^n- V<j)±(x) = T^r(x ) + K [r)(x ) 
lim ± n x V(j)±(x) = Kx[r](xo), 

K [r](x ) = [ |^(«o - x)r(x)dA(x) 

Ki[r\(x ) = n x V J g k (x - x) r(x) dA(x) , 

Kq is an integral operator of order — 1 and K\ is an integral operator of order 0, 
which is defined in a principal value sense. 
For the vector potential we have 

lim n ■ A±{x) = K 2 , n \j]{x Q ), 



(29) 



lim n x A±(x) = K 2) t[r](x ) 



where 



K2,n[j]( X o) 



g k (x - x) (n • j(x)) dA(x) , 
g k (x - x) (n x j(x)) dA{x) 



i^2,n an d Ki t are both integral operators of order — 1. 
The vector potential also satisfies 



where 

K 3 [j](x ) 



lim n • V x A±(x) = K 3 [j](x ) 
lim n x V x A±(x) = ±\j{x Q ) + K 4 [j](x ) 



J Vg k (x - x) ■ (j(x) x n ) dA(x) , 



(30) 



Vg k (x - x) (j(x) ■ n ) 



dgk 
dn n 



x - x)j(x) 



dA(x) . 



K3 is an integral operator of order 0. K4 is an integral operator of order — 1. 
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Proof. These results follow from classical potential theory, the observation that the 
kernels in K^ n and are weakly singular, and the fact that, at the singular point 
xq = x in K4, no is orthogonal to j(x)). □ 

Remark 5. We have abused notation slightly in the preceding Lemma. The opera- 
tors are functions of the Helmholtz parameter k. When the explicit dependence is 
relevant, we will occasionally write K (k), K\{k), K±{k) instead. 

The limits of the anti-potentials A m and 4> m are analogous. Recall, however, 
that we have chosen not to work with j and m as unknowns, but rather the gener- 
alized Debye sources r, q complemented by the harmonic vector fields, j and m 
are computed from 

j = Vrf + V r x (ntf m ) + j H 
m = n x j 

where * and ^ m satisfy the Laplace-Beltrami equations (|25T > with r, q viewed as 
source data. 

Lemma 4. The integral operators K2, K3 and K4 are all of order —1 or —2, and 
hence compact, when viewed as operators acting on r, q. 

Proof. This follows immediately from Lemma|3]and the fact that j, m are of order 
— 1 in terms of (r, q). □ 

From Lemma|3l we obtain the following jump relations: 

Corollary 2. Suppose that the fields E,H, are defined in terms of potentials and 
anti-potentials. Then they satisfy 

n-{E + -E_) = r n x (E+ - E_) = -m 

n ■ (H + - H ) =q nx (H + - H ) = j ' { ' 

4 The Maxwell Equations with Normal Components Spec- 
ified 

We note that, for k G C+, the fundamental solution 

e ik\x\ 

9k{x) = — j— 1 (32) 

is outgoing; that is, it satisfies the Silver-Miiller radiation condition. Thus, all of 
the corresponding potentials defined over bounded regions are outgoing as well. 
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Theorem 4. Let E and H be outgoing fields represented in terms ofDebye sources 
(r, q) and currents j , rn. Then the limiting values of their normal components are 
given by the following integral representation for Xq £ T. 

±\I-K Q \(r\ (ikK 2>n -K 3 \fj\ 
t\I + K ) \qj ^\ K 3 ikK 2 J\m)> 

where I denotes the identity operator. If we assume that j = j R (r,q,k) and 
m = n x j , then these are Fredholm integral operators of the second kind in the 
generalized Debye sources (r, q). As k tends to zero, this representation converges 
to 

j5t(s::)-( ±i/ -. J6W #:*m)Q - 

The equation [—1 + 2Ko(0)]f = h is uniquely solvable for all h, and the equation 
[I + 2X0(0)]/ = hfor all h of mean zero. 

Proof. The equations follow from the formulae in the previous section. The solv- 
ability properties of / ± 2Kq(0) are classical and can be found in [H. □ 

Definition 4. We let AA ± (fc) denote the operator on the right hand side of (l33l >. 
with J\fg(k) the first row and M%(k) the second. 

If we seek to impose the boundary conditions E n |p= /, H ■ n |p= h, then 
we obtain the following system of equations, which is analytic in k. 



(±\I -K \ fr\ (ikK 2:n -K 3 \fj\ ff\ 

V Tp + KoJ \q) + V K 3 ikK 2tn J\m) \h) ^ 

When j = jji(r,q,k) and m = n x j are obtained from the Debye sources via 
(|25T >. we denote the left-hand side of (I35T ) by J\f ± (k) (r, q) and observe that it is a 
compact perturbation of the operator J±(0) where 

Theorem 5. Let (/, h) E A^r.o- For k £ E + , a discrete set in the complex plane, 
the equation 

Ar + (k)(r,q)=(f,h) (37) 



15 



has a unique solution. The outgoing solution of the THME(k) defined by this data 
satisfies 

E+-n\ r =f H + n \ r = h. (38) 

Proof. We know from Theorem @] that the operator J+(0) is invertible. By analytic 
Fredholm theory, therefore, we know that there is a discrete set E + 6 C so that 
for k £ E + , N + (k) is also invertible. The result now follows from the fact that 
(r, q) £ Mr,o (LemmaQ} and that (/, h) £ Mr,o (Lemma©, so that 

M + (k) : M r ,o -> Mr,o- 

□ 

Corollary 3. When T is simply connected, the Fredholm equation < \37i provides a 
unique solution to the THME(k)for any k £ C+. 

Proof. This follows immediately from Theorem [5] and the fact that E + lies in the 
complex lower half -plane, as follows from Corollary [4] in Section [7] □ 

Remark 6. The problem of solving the THME with specified normal components 
was previously analyzed by Gulzow lfl2l . who constructed a hypersingular inte- 
gral equation method. Using generalized Debye sources instead leads to a well- 
conditioned integral equation of the second kind. In the multiply connected case, 
Theorem [5] shows that the problem of scattering with normal components speci- 
fied is invertible, if the solution is sought with zero projection onto the harmonic 
vector fields. The extra conditions in Kress' result (Corollary [U force uniqueness 
on j H by specifying g conditions on the circulation of E and g conditions on the 
circulation of H, where g is the genus of T. In fact, any set of 2g conditions that 
uniquely determines the harmonic component of the current, j H , will suffice, (see 
Theorems [9] and fTTT) . 

5 The Perfect Conductor 

We turn now to the problem of scattering from a perfect conductor, which requires 
an analysis of the tangential components of E and H. These are easily expressed 
in terms of potentials using Lemma[3] 
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Theorem 6. Let Xq 6 V. The limiting values of the tangential components of 
E, H, in terms of potentials and antipotentials, are given by 



X 




(r\ 








ikK 2 , t 



-K A 
ikK 2 ± 



) 



q 

3 



(39) 



where K\, Ki t, K4 are defined in Lemma\3\ 

Definition 5. In the sequel we denote the operator on the right hand side of (l39l as 
T ± (k). We use T^(k) to denote the first row, and T^{k) to denote the second. 

For scattering from a perfect conductor, let us recall that both ([8]) and Q must 
hold on T. As noted in the introduction, the EFIE approach involves imposing ([8]) 
using only the classical vector and scalar potentials, A and <fi, with the physical cur- 
rent J as the unknown. This leads to an integral equation on T with a hypersingular 
kernel that has interior resonances and suffers from low-frequency breakdown. To 
avoid these difficulties, we turn again to the Debye sources. A nonstandard feature 
of our approach is that we will extract only one scalar equation from the tangential 
conditions on the E field and couple it to the normal condition © satisfied by H . 

5.1 The hybrid system 

The operator defining the tangential components of E± is given by 



If we restrict to j = jji(r, q, k), and m = nxj , then, acting on (r, q) G .Mr,0) the 
only term of non-negative order is — K\r. We use Tg(k)(r, q) to denote this oper- 
ator restricted to this subspace of data. In order to recast K\ in (l40l ) as a Fredholm 
operator of the second kind, it is convenient to multiply it by a left parametrix, 
based on the following standard result. 

Lemma 5. Let T C ffi 3 be a smooth, connected closed surface, let xq £ T, and 
let Go denote the single layer potential operator based on the Green's function for 





Q 
3 



(40) 
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the Laplace equation: 

G Q [f](x Q )= J 9o (x-x )f(x)dA(x), (41) 
and let (ft denote the usual scalar potential with density r. Then 

lim G [A T (ft](x) = jr(x Q ) + K^xo) (42) 

where K\ is a compact operator. 

Proof. To see this, we observe that the surface Laplacian satisfies the identity 



A r (ft = Acft- 2H-- 



dn dn 2 ' 

where H denotes mean curvature (see, for example, (211). Since Acft = -k 2 (ft (by 
construction), we have 

n deft d 2 (ft 

A r = -k 2 <ft - 2Hjf - — T . (43) 
on on z 

The composition of Go with the first two terms on the right-hand side of (l43l are of 
order —2 and —1, respectively, hence compact. It is a classical result (a Calderon 
relation) that the composition of a single layer potential with the second normal 
derivative of (ft yields a compact perturbation of \l II2T1 : 

Jh^ ± G o [0](z) = \r{x Q ) + D 2 [r](x ) 



where 



D[r](x ) = J ^-(x -x)r(x)dA(x). 



Here n T denotes the normal at x, and -?p— denotes the normal derivative at x. D is 
the usual double layer potential, which is the adjoint of Kq. Since D is an operator 
of order —1, the result follows. □ 

Lemma 6. Let T C M 3 be a smooth, connected closed surface, let xq £ T, let Gq 
denote the single layer potential operator (HZ]). Then 



lim^ G Vr • [n x T±(k)] ( ' n ) = -^(xo) + iV x (fe) ^ ) . (44) 



where Ni(k) is an analytic family of operators of order -I. 



18 



Proof. The result follows from Definition [51 the fact that 

V r • n x Ki[r] = -A V (f>[r] , 
and the preceding Lemma. □ 
Recalling that 

j R (r,q,0) = 0, 

we see that 

G V r • [n x T±] (0) Q = G [A r G r] . (45) 

If T has M components, then the nullspace of this operator is M -dimensional. It 
is generated by functions r such that G$r is constant on each component of F. As 
a consequence of Theorem 5.7 in (H, it follows that this nullspace only intersects 
Mr,o at 0. 

Definition 6. Taking the integral operator in Lemma [6] and the integral operator 
Nh from DefinitionHJ we define Q ± (k) as the following hybrid system of integral 
operators: 

As divergences, the range of Vr • [n x T^] (k) consists of functions of mean zero. 

Proposition 1. The family of operators Q ± (/c) is analytic in k and Fredholm of the 
second kind. There is a discrete subset F + C C, such that Q + (k) : Mr,o —> Mr,i 
is invertible for k ^ F+, where Mr,i is the l?-closure of 

{(G r,q) : (r,q) € M r ,o}- (47) 

Proof. The analyticity statement is immediate from the formula. Examining Q 1 * 1 (k) 
we see that 

where is an analytic family of operators of order —1. Thus, Q ± (k) is Fredholm 
of the second kind. The last statement follows from the facts that Q + (0) : Air,o — * 
Mr,i is invertible and Q + (k)Mr,o C Mr,i- □ 

We may now state our principal result in the simply connected case. 
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Theorem l.IfY is simply connected, then F + is disjoint from the closed upper 
half plane. Thus, the integral equation 

e + W ( r 8 ) = ({) (49) 

provides a unique solution to the scattering problem from a perfect conductor for 
any k in the closed upper half plane. Here, 

f = G V r -[nxnx E in ], h = n-H in \ Tr . (50) 

This is proved as Theorem [12] in Section 1770 We leave the discussion of applying 
our method in the non-simply connected case to Section 17711 

5.2 Low Frequency Behavior in the Simply Connected Case 

The representation of solutions to the THME(&;), using data from Mr,o © W 1 (T), 
afforded by (l83l behaves well as the frequency tends to zero. In the simply con- 
nected case we only have data from A^r.o- As k tends to zero, this space of so- 
lutions tends to the orthogonal complement of the harmonic Dirichlet fields, that 
is outgoing harmonic fields, with vanishing tangential components on 60, so that 
E and H are recovered from — V</> and —VcpM alone. This is proven, along with 
the multiply connected case in Section 17721 It is therefore apparent that Q + {k) 
provides a means for finding and representing solutions to the perfect conductor 
problem, which has neither interior resonances, nor suffers from low frequency 
breakdown. 



6 The Exterior Form Representation 

For the remainder of this paper we represent the electric field E as a 1-form and 
the magnetic field H as a 2-form, rj. This choice is explained in Appendix [A] If 

E = eid Xl + e 2 d X2 + 63^3 and H = h\d Xl + h 2 d X2 + h 3 d X3 , (51) 

then 

£ = &\dx\ + e 2 dx2 + e^dxs and 

rj = h\dx2 A dx% + h^dx^ A dx\ + h^dx\ A dx 2 - 

If (-, •) denotes the Euclidean inner product, i.e. the metric on M 3 , then £ is defined 
by the condition that, for every vector field V we have: 

(V,E) = £{V). (53) 
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That is, £ is the metric dual of E, and (with ★ the Hodge star-operator, see 
Remark[8]and Appendix |A.2b is the metric dual of H. 

The curl-part of time harmonic Maxwell's equations takes the form: 

d£ = ikr) d*rj = -ik£. (54) 

For k ^ 0, these equations imply the divergence equations, which take the form: 

d*£ = dr) = 0. (55) 

An outgoing solution to the Helmholtz equation on 1 -forms satisfies the analog of 
the Silver-Muller radiation conditions: 



i&d$, ~ d*£ x ■ dx — ik^ = o ( - — - 



(56) 



where x = -p^. A magnetic field, rj is outgoing if *r] satisfies (I56I ). 

The standard integration by parts formula in electromagnetic theory is derived 
by considering the L 2 -norm of the quantity in the radiation condition. We let Dr 
denote the ball of radius R, centered at 0, Sr = bDji, and CIr = Dr n £1. Since 
the quantity in (l56l ) is o(|x| _1 ) it follows easily that 



lim 

R— >oo 



\i & d£ - d*£ x ■ dx - ik£\\ 2 dA = 0. 



(57) 



We expand the integrand to obtain that 

\ixd£ -d*£x- dx\\ 2 + 



lim 

R— too 



dA+ 



( 



-2 Re 



k J (H,i A d£)-(i A td*OdA 



\ Sr 



(58) 



Using Green's formula, we can replace the second integral with 



2Im(/c) 



+ IK£f + l^ril^f I dV+ 



21m 



k J [(&i n dt)-(i n &d*g)] dA 



(59) 



Here we use n to denote the inward pointing normal field along F. 

Combining (I58T ) and (|59l ) we obtain the standard integration by parts formula 
for outgoing solutions to the vector Helmholtz equation: 
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Lemma 7. If £ is a 1-form defined in Q, satisfying A£ + k 2 £ = and the radiation 
condition (1561 ), with \m{k) > 0, then 

lim [2Im(Jfe) / [\\d£.\\ 2 + \\d*£\\ 2 + \k\ 2 \\£\\ 2 ]dV+ 

R^oo \ J 
V Mr 

j [IM£ - <rtx- dx\\ 2 + \k\ 2 u\\ 2 ]dA^j = (60) 

r 

Remark 7. If (£, 77) is a solution to the THME(A;) in fi, then the outgoing radiation 
condition can be rewritten as 



1 



(61) 



in agreement with ©. 



6.1 Uniqueness for Maxwell's Equations 

If (£, 77) is a solution to the Maxwell system, then d*£ = and the boundary term 
in (l60l ) reduces to 

-21m fife J {£,ind$)dA^ . (62) 

Let v denote a 1-form defined along T, which restricts to zero on TV and is nor- 
malized by v(n) = 1. 

Remark 8. We use * to denote a Hodge star-operator, see Appendix IA.2I In much 
of the paper we need to distinguish between the Hodge star-operator acting on 
forms defined on M. 3 and that acting on forms defined on surfaces in M 3 . We denote 
the M 3 -operator by *3, and a surface operator by *2- Which surface is intended 
should be clear from the context. 

If T is a smooth closed surface in M 3 , which bounds a region D, then it ob- 
tains an orientation from its embedding into R 3 : let n be the outward pointing 
unit normal vector, and X\, X2 a local oriented orthonormal frame for TV. We let 
ui,U2, v, be the local co-frame for T*IR 3 |r, dual to X\, X2, n. Note, in partic- 
ular that, uj\(n) = u>2(n) = 0. We say that the frame (Xi, X2,n) (or co-frame 
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(ui,u)2, v)) is adapted to V. The 1-form that is the metric dual of the vector field 
aX\ + bX2 is auj\ + bu)2 ■ 

The volume form on M 3 and area form on T are given locally by 

dV = uj\ A L02 A v dA = i n dV = u\ A uj 2 . (63) 

In terms of the adapted frame, the Hodge star-operator on M 3 is 

★3I = uj\ A u>2 A v *3 wi A UJ2 A ^ = 1 

= (x>2 A V *3W2= — Wl A V *3 V = LJi A Li>2 (64) 
*30>i Al^ = — U>2 *3 W2 A V = UJ\ *3 Wi A W2 = 

The Hodge star-operator on the surface (oriented as the boundary of D) is given by 

★2I = k>i A a>2 *2 A c<J2 = 1 

(65) 

★2W1 = W 2 *2 W 2 = — Wl. 

It is useful to note that if q is a 1-form defined on F, then *|q = — q. If v is a vector 
field tangent to V and u, its metric dual, then *2u; is the metric dual ofnxi). 

To emphasize the distinction between an exterior form acting on TR 3 restricted 
to a surface rcl 3 and the restriction of this form to directions tangent to T, we 
sometimes use a |r to denote the former notion of restriction, and a \tt, the 
latter. For a 1-form, a, represented along T in the adapted co-frame by a |r= 
au\ + buj2 + cv, we have 

a \tf^ awi + bu)2- (66) 

We denote this latter restriction by a t - We use the notation d? to denote the exterior 
differential acting on forms on T. 

For a 2-form (3 \r= 00*1 A v + buj2 A ^ + clji A U2, we have 

/3 |Yr^ cwi A c^- (67) 

For a 3-form 7, dimensional considerations imply that 

7 trr= 0. (68) 

If / is a 0-form, or scalar function, then 

/ Trr= / \r . (69) 
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The definition of the inner product on forms, the fact that (i>, i n d$) = 0, and 
the definition of *2 on T, imply the identity 



2ImlkJ (£,i n d£)dA J = -21m ( k J \ r ] A* 2 M£ 




(70) 



Using this identity and Lemma [7] we can prove the basic uniqueness theorems. 
Note that £ is a 1-form so £ t = £ \tt corresponds to the tangential components 
of E, and i n £ |r, the normal component. The magnetic field rj is a 2-form and 
therefore r\ \tf gives the normal component, and [i n f]]t = inV |tt, gives the data 
in the tangential components, corresponding to n x H. 

We restate the classical result (TheoremQ]) that an outgoing solution to THME(A;) 
is determined by either the tangential components of the electric or magnetic fields 
in the language of forms: 

Theorem 8. Suppose that (£, rj) is an outgoing solution to the THME(k), in Q,for 
a k ^ 0, in C+. If either £ t , or [i n r]]t vanish, then the solution is identically zero 
in £1. 

Kress' result (Theorem© on the normal components of (£, rj) is restated as 

Theorem 9. Suppose that (£, rj) is an outgoing solution to the THME(k), in Vt,for 
ak ^ 0, in C+. If every component of the boundary ofVt is simply connected, then 
the solution is determined by the normal components i n £ fr ond rj |Yr • If the 
sum of the genera of the components ofT equals g > 0, then there is a subspace of 
outgoing solutions to THME(k) with 

i n £ [r= and r\ \ Tr = (71) 

of dimension 2g. 

Remark 9. As noted above, we call solutions to THME(fc) that satisfy dTTT ) k- 
Neumann fields, and denote the space of such solutions by Hfe(fi). Here we give 
a bound on dimHfc(T2), and a novel description of the additional data needed to 
specify the projection into this space. Later in the paper we give a new proof that 
dimTik{Q) = 2g, for k in the closed upper half plane, C+. In this regard, the case 
k = is classical. As noted above, this result was proved in 031. 

Proof. Suppose that i n £ and rj |r both vanish. Let a = £ |r and (3 = ^377 |> • 
The usual properties of the exterior derivative, the hypothesis 77 fr= 0, and the 
equation d£ = ikrj imply that 

d v a. = 0. (72) 
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We can rewrite d*rj = —ik^ as d*3 r\ = —ik *3 The hypothesis, i n £ = now 
implies that 

d r /3 = 0. (73) 
A calculation using a co-frame adapted to T shows that 

*2j3 = -M]rr, (74) 
see (l97l) . The equation d£ = ikr\ implies that 

\i n d£\ \r=-ik* 2 p. (75) 
We can therefore express the right hand side of d70l) as 



2\k\ 2 Re / a A/3 . (76) 



If T is simply connected then the equation dpc* = implies that a = dru. As 
dr(3 = as well, a simple application of Stokes formula shows that 



/ 



d TU A (3 = 0. (77) 



This completes the proof of the theorem when T is simply connected. 

For the general case, let #j R (r) denote the De Rham cohomology group with 
dim Hl R (T) = 2g ^ 0. We show that the space of solutions with vanishing normal 
components, for which the integral in (|76"1) is non-vanishing, depends on at most 2g 
parameters. Using the wedge product, we define a pairing, W, on closed 1 -forms: 

W{r],u) = J i] Auj. (78) 
r 

If dri] = and uj = dru, then, as noted above, Stokes theorem implies that 

W(rj,u) = 0. (79) 

Hence W defines a skew-symmetric form on if3| R (r), which is well known to be 
non-degenerate. As the dim ^ R (T) = 2g, this observation completes the proof of 
the fact that dimWfc(O) < 2g. If the image of either a or (5 in _£fj R (r) vanishes, 
then W(a, (3) = 0, which implies, as above, that the solution in 17 is zero. □ 
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From this proof we see that the image of either a or (3 in H 1 (F) provides data spec- 
ifying a £>Neumann field. This should be contrasted with the data used by Kress, 
given in equation (|2"7T ). As dimWfc(S7) = dim if 1 (r), the maps from to 
i7 1 (r) defined by (£,f]) i— > a and (£,rj) i— > /3 are both isomorphisms, when 

We complete this section by proving Lemma |2j which shows that the normal 
components of (£,rj), a solution to THME(fc) with k ^ 0, have vanishing mean 
value over every component of the boundary. While superficially this might appear 
analogous to the fact that the normal derivative of a harmonic function in a bounded 
domain has mean value over the boundary, it is actually an elementary consequence 
of the equations themselves and Stokes' theorem on a closed surface. 

If k ^ 0, then the Maxwell equations, (l54l ). imply that 

*3 £ tr= ^ d r[*3V \t] V \r= ^d r [£ \r\. (80) 

As 

*s £ \r= i n £dA (81) 

the relations in (l80l imply that these forms are exact and therefore Stokes' theorem 
implies that a solution of THME(fc), with k ^ satisfies: 

j i n £dA = j V = (82) 

for m = 1, . . . , M. Note that this is true whether the limit is taken from Q or D. 
This completes the proof of Lemma|2j restated as 

Proposition 2. Let (£, rj) be a solution to the THME(k)for a k ^ 0, in a region 
Gel 3 , with a smooth bounded boundary. The normal components (i n £, i n *3 rj) 
have mean zero over every component ofbG. 

6.2 Potentials and Boundary Integral Equations 

We now re-express (fl~6l ) in terms of exterior forms. Assuming, as before, that 
the time dependence is e~ twt , the permittivity is e, the permeabilty is /i, and k = 

u-y/e/J; we set 

£ = (ika — d<p — d*a m ) rj = (da + ika m + d*$ m ), (83) 

where <f> is a scalar function, a a one form, a m a two form, and <& m = 4> m dV, 
a three form. This representation is quite similar to what one obtains using the 
fundamental solution for the Dirac operator d + d*, see 0. 
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In order for (£, t}) to satisfy the THME(fc), the potentials must satisfy: 

d*a = —ikcp da m = ik& m . (84) 

As before g k {x — y) denotes the outgoing fundamental solution for the scalar 
Helmholtz equation, with frequency k. As discussed in the introduction, all of the 
potentials are expressed in terms of a pair of 1 -forms j, m defined on T, though 
in the end, we do not use j and m as the "fundamental" parameters. When we 
express these 1-forms in terms of the ambient basis from M 3 , e.g., 

3 = h{x)dxi + j 2 (x)dx 2 + j 3 (x)dx 3 , (85) 

we normalize with the requirement 

i n j=j(n) = 0. (86) 

These 1-forms are the metric duals of the vector fields, tangent to T,' previously 
denoted by j and m. 

The "vector" potentials are given in terms of surface integrals by setting 



ol = I 9k{x - y)[ji(y)dxi + j 2 {y)dx 2 + j 3 (y)dx 3 ]dA(y) 

(87) 



r 



CX-m = *3 



9k{x - y)[mi(y)dx\ + m 2 (y)dx 2 + m 3 {y)dx 3 ]dA(y) 



x 



Using the equations in d84l ) we obtain the form of the potentials defining and 
$m = 4> m dV, letting 



4>(x) = I 9k(x - y)r(y)dA{y) 

(88) 



r 



4>m{x) = I g k {x -y)q{y)dA(y), 



where we let 



-r^r *2 j = rdA \dr * 2 m = qdA. (89) 

ik ik 

The scalar functions, (r, q) are, as before, the Debye sources. From this defini- 
tion, and Stokes' theorem we see that the mean values of r and q vanish on every 
connected component, Fj, of T, 

rdA = / qdA = (90) 
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This proves Lemma [T] It is necessary for the conditions in d89l to hold, and thus, 
for (£, rj) to satisfy the Maxwell equations. 

As before, we let Air,o denote pairs of functions (r, q) defined on T with mean 
zero on every component of T. If we assume that m = *2j, (as we usually do), 
then, taking account of the fact that the (negative) Laplace operator on 1 -forms is 
given by — Ai = d^dr + a?r<if , we obtain the relation: 

Aij = ik[d r r - * 2 d r q]. (91) 

If the components of T are all of genus zero, then equation d9TT ) always has a unique 
solution. If T has positive genus components, then one needs to deal with the null 
space of Ai. 

If Hl R (T) / 0, then the nullspace of Ai, H 1 ^), which agrees with the space 
of solutions to 

d r a = d^ot = (92) 

is isomorphic to iJl R (r). These are the harmonic 1-forms. The right hand side 
in (I9TI ) is orthogonal to W 1 (r), and hence lies in the range of Ai. Let i?i denote 
the partial inverse of the Laplacian on 1-forms, with range orthogonal to W 1 (r), 
and set 

j R (r, q, k) = ikRi [d T r - * 2 d r q] , (93) 

Because the ranges of dr and d^ are orthogonal to the null space of Ai, this equa- 
tion is solvable whether or not r and q satisfy the mean value condition. Note that 
the solution to (l93l) tends to zero as k — > 0. 

Using the relations Aidr = dr&o, and A2 = *2Ao*2> we can re-express j R 
in the form: 

j R (r, q, k) = ik[d r Ror - * 2 drRoq\- (94) 

Here Ro is the partial inverse of Ao, which annihilates functions constant on each 
component of V and has range equal to the set of functions with mean zero on each 
component of T. Equation d94l shows that this approach to representing solutions 
of Maxwell's equations in terms of the pair (r, q), only requires an inverse for the 
scalar Laplacian on T. 

For any (r, q), the solution space to (I9TI ) is isomorphic to W 1 (r). Adding a 
harmonic 1-form to j R does not change r and q, though it changes the fields £ and 
rj, and plays a central role in the discussion of Hfe(fi). If g ^ 0, then the space of 
outgoing solutions to THME(A;) is parameterized by Air,o © T~l l (T). So given data 
(r,q,j H ) we often speak of the solution to the THME(/c) "defined" by this data. 
If j H is missing, then it should be understood to be zero, i.e. the solution "defined 
by (r, q)," is the solution defined by (r, q, 0) in the sense above. 
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6.3 Boundary Equations and Jump Relations 



As with the vector field representation, we can take the limits of £ and rj in (I83T ) as 
the point of evaluation approaches T, and obtain boundary integral equations for 
the normal and tangential components of these forms. Indeed there is no necessity 
to rewrite these equations, we simply use 7^ (k) (r, q,j, m) , (k) (r, q, j , m) to 
denote the limiting tangential components and the limiting normal components are 

(k) (r, q, j , m) , M^(k) (r, q, j , m) . Keep in mind that, in the vector field rep- 
resentation, the tangential components are represented as the limits of n x £ and 
n x H, which correspond to * 2 £ t and *2([*3*?]t)> respectively. For consistency, 
we use , to denote the boundary values of these quantities. As before, + in- 
dicates the limit taken from Q and — the limit taken from D. If j = j R (r, q, k) and 
m = *2j\ then we omit them from the argument list, e.g., we use the abbreviated 
notation T±(k)(r, q) to denote T±(k)(r, q,j R (r, q, k),* 2 j R (r, q, k)), etc. 

Below we use the jump relations to prove a uniqueness theorem. So it is useful 
to reformulate them in the form language. The magnetic field is represented by 
the 2-form, 77 = h\dx 2 A dx% + h 2 dx^ A dx\ + h^dx\ A dx 2 , so that ^377 = 
hidxi + h 2 dx 2 + hsdxs. The most direct way to define the normal and tangential 
components of r/ along T is as i n *3 r]± and (*3t]±)t- In this formulation the jump 
relations then take the form 

i n (£+ -£-) = r (£+ - = * 2 m ^ 
i n (*3V+ ~ *3»7-) = Q (*3*?+ - *3V-)t = -*2j 

It is also useful to calculate the relationship between i n rj and (*3rj) t . In an adapted 
frame (loi,lo 2 ,u), We have 



77 = au>i A uj 2 + bu A U\ + cv A uj 2 (96) 



and therefore 



I i-nXI = bu>i + COJ 2 

<. , , which implies that * 2 [i n V\ = {*3V)t- (97) 

[(*3»7)t =buj 2 -cwi, 



7 Uniqueness for the Tangential Equations 

Suppose that there is a k 6 C + \ {0}, and non-trivial data (r, q) G .Mr,o an d j, 
satisfying d89l , with m = * 2 j, so that 

T+(k)(r,q,j,m) = 0. (98) 

Let (£ ± , T7 ± ) be the solutions to the Maxwell equations defined by this data in the 
complement of T. By their definition it is clear that the tangential components of 
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(101) 



£ + vanish along T. As Im(fc) > 0, the solution in Q is outgoing and Theorem [8] 
shows that (£ , , rj , ) = (0, 0). The jump relations, (1951 ), and the fact that m = *2jj 
allow us to determine the tangential boundary data for (£_,rj_) : 

£_ \tt=3 inV- \tf = 3 ■ (99) 

It is not difficult to see that the boundary condition on the Maxwell system in D, 
implied by these relations, 

$- \tt= inV- \rr, (100) 

is not formally self adjoint! 

Observe that d*d£_ = k 2 £_ and i n d£_ |r= iki n rj_ |r • Using a standard 
integration by parts formula, we obtain: 

/ \d£_,dt_)dV = j(d*dd_,£_)dV + j \i n d£_,£_)dA 

D D bD 

= k 2 j \t_,Z_)dV + ik J(i nV _ \ r ,t_ \r)dA. 

D bD 

Combining this with d99l > gives: 

-ik j (j,j)dA = k 2 j (£_ t £_)dV- f(dZ_,dt_)dV. (102) 
bD D D 

We can rewrite this relation as 

- aik = bk 2 - c, (103) 

where a, b, c are non-negative real numbers. If b or c vanishes, then it is clear that 
£_ = 0. If a = 0, then j = 0. For a countable set of real numbers {kj}, there exist 
non-trivial solutions to the equations 

d*d£_ = k 2 £_ d*£_ = £ \ TbD = 0. (104) 

In the present circumstance, however, (r, q) are generalized Debye sources and 
therefore ^ 

rdA = —dr *2 3 and qdA = —drj- (105) 

ik ik 

If a = 0, then all the boundary potentials vanish, and therefore £ = as well. 
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Using the quadratic formula we see that 



— ia ± y/Abc — a 2 
k ± = - . (106) 

As a,b,c are all positive, (11061 ) shows that Imk± < 0. This argument applies, 
mutatis mutandis to T T j~(k). Formula (11061 ) and the discussion above complete the 
proof of the following theorem: 

Theorem 10. Assuming that m = *2j> an d (r, q) satisfy d89l ), then, for Im k > 0, 
k 0, the nullspaces of both T^~(k) are T^(k) are trivial. 

In the case that V is simply connected this implies that E + , the exceptional set 
for Af + (k), is disjoint from the closed upper half plane. 

Corollary 4. If every component ofT is simply connected, then, for k with Im k > 
0, the Fredholm operator M + {k) is an isomorphism from TWr.o to itself. For such 
k, the rows ofT + (k) are also surjective and hence isomorphisms. 

Proof. If T is simply connected, then any 1-form j on T has a unique representa- 
tion as j = jn(r, q, k). We can therefore regard N + {k) as a Fredholm system of 
second kind for the normal components of rj + ) in terms of (r, q). In this case, 
Theorem [9] implies that a solution (£ + ,rj + ) of THME(fc), with vanishing normal 
components is identically zero in f2. Hence (£_, r)_) satisfy (|99l ), and we can there- 
fore apply the argument leading up to Theorem[l0]to prove that E + is disjoint from 
the closed upper half plane. The Fredholm alternative then implies that J\f + {k) is 
also surjective. The surjectivity of the rows of T + {k) is now immediate. □ 

Remark 10. The poles of the scattering operator for the Maxwell system, defined 
by a self adjoint boundary condition on T, lie in the lower half plane. Nonetheless, 
it appears that the eigenvalues for the non-self adjoint boundary value problem 
defined by (1 1001 ) are unrelated to these poles, but are simply interior resonances, 
familiar from more traditional representations of solutions to Maxwell's equations 
(see the Introduction and Remark ITSb. The non-self adjointness of this B VP places 
the interior resonances in the lower, non-physical, half plane. This leads, in the sim- 
ply connected case, to numerically effective algorithms for solving the THME(/c), 
which do not suffer from the instabilities caused by interior resonances in the phys- 
ical half plane. 

In the non-simply connected case we have the following theorem assuring the 
existence of A;-Neumann fields. 

Theorem 11. For k € C+, the space of k-Neumann fields has dimension exactly 
2g, and the rows ofT + {k) are surjective. 
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Proof. For k G C+, the solutions to THME(fc), defined by 

C H = {(0,0,3 h,*23h) ■ iffe^(r)} 
have a trivial intersection with those defined by data in 

Cr = {(r,q,j R (r,q,k),* 2 j R (r,q,k)) : (r,q) E Mr,o}- 

The solutions defined by data in Ch may not themselves be fc-Neumann fields. 
To find solutions in Hfc(f2), we first use an element of Ch to construct a solution 
(£o+> ? ?o+)- If ^ ^ then we can solve 

A/"+(A0( r> ) = (> e °+ V (107) 

and denote the solution of the THME(£;) defined by this data in by (£ 1+ ,rj 1+ ). 
By Theorem [TOJ the difference 

(£jv+, Vn+) = (€o+,Vo+) ~ (£i+,Vi+) ( 108 ) 

is a non-trivial A;-Neumann field. These solutions depend analytically on k G <C + \ 
E+. As (£jv+> Vn+) i s a non-zero solution to the THME(A;) with vanishing normal 
components, Theorem [9] shows that the cohomology class of £]y +t must be non- 
trivial. Thus for k ^ E + , dim Wfc(f2) is at least 2g. On the other hand, the proof of 
Theorem [9] gives the upper bound dim Ti^VL) < 2g. Proving the theorem in this 
case. 

Now suppose that kj G E + n C+. This means that there is a non-trivial, fi- 
nite dimensional space of data V^. C A^r,o> which defines ker J\f + (kj)(r, q). 
Let (£±,»7±) denote the solution of the THME(A;) defined by a non-zero pair 
(r, q) G VL ■ The restriction of £ + to TT defines a cohomology class in H^ R (T). 
If this class is trivial, then Theorem [9] implies that the pair (£, ,tj.) are identi- 
cally zero. In this case, Theorem [lOl implies that Imkj < 0, contradicting the 
assumption that it lies in C+. This establishes that each non-trivial pair in 
defines a non-trivial /c, -Neumann field, thus a subspace of of dimension 

d = dim Vkj . 

The Fredholm alternative implies that the equations for the normal compo- 
nents: M + (kj)(r,q) = (f,g) are solvable for pairs (/, g) G Mr,o satisfying 
exactly d linear conditions. This means that within Ch there is a subspace of di- 
mension at least 2g — d, for which the normal components can be removed, as 
above. We therefore get another subspace, c Hfe .(fi), of dimension at least 
2g — d. As V] ij has a trivial intersection with the data defining U^., Theorem [lOl 
implies that these two subspaces of Ttk have a trivial intersection. The lower 
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bound on the dimension of £7j,. and the upper bound on dim 7^. (fi) imply that 
dim Ukj = 2g — d; this completes the proof that dim 7^(0) = 2g, for all k € C+. 

For k ^ Theorem [9] combined with the fact that dim 7^(0) = 2g shows 
that the rows of T + {k) are surjective. If k G E + n <C + , then the range of M + {k) 
has codimension exactly d = dim . On the other hand, there is a d-dimensional 
space of data in TL l (T) for which the normal components span a complement to 
that in \mj\f + (k). Once again we can find an outgoing solution to the THME(fc) 
with any specified normal components. Combined with the fact dim 7^(0) = 2g, 
Theorem [9] again shows that the rows of T + {k) are surjective. □ 

In the course of this argument we established: 

Corollary 5. For k € C + \ {0}, the map from Hfe(O) to H^iT) defined by 

is an isomorphism. 

Remark 11. It was shown by Picard that for each k with non-negative real part 
there are families of interior ^-Neumann fields, that is non-trivial solutions to the 
THME(A:) in D with vanishing normal components, see Il23~ll24l . For most values 
of k there is a 2g-dimensional family. In lfI31 Kress showed that there is a countable 
set of positive real numbers {kj}, with kj — ► oo, for which there are non- trivial 
interior A; -Neumann fields with vanishing circulations. 

This theorem really asks more questions than it answers: 

1. If A; = 0, then Hodge theory essentially implies the existence of the 0- 
Neumann fields. For k ^ 0, what is the reason for the existence of the 
A;-Neumann fields? A possible explanation might go along the following 
lines: In Appendix I A.4[ we express the THME(A;) in the form: 

LkiA + ri) = (d + d*-ikA)(£ + rj) = 0. (109) 

Suppose that acting on divergence free, outgoing data, which satisfies 
in£ = V \r= is in some sense a Fredholm family. The boundary condi- 
tions defining the formal adjoint, are £ t = 0, (i n v)t = 0. Theorem [8] 
implies that the nullspace of [L^]* is trivial for k G C+. The nullspace of 
Lfc is exactly 7Yfc(f2). If these operators are a Fredholm family, then the con- 
stancy of the Fredholm index would imply that 

Ind(L fc ) = dimWfe(O) = 2g. (110) 

It is not obvious, however, on what space the range of is closed. 
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2. Does E + have a non-trivial intersection with Im k > 0? If so, what is the 
physical significance of these numbers? As noted above, Kress proved that 
there is a countable set of positive real numbers for which there exist interior 
fc-Neumann fields, with vanishing normal components and circulations. Are 
these numbers in any way related to E + l 



7.1 The hybrid system using forms 

The operators defining tangential component of are given by 

(r\ 



r±(fc) 



THt + (- Kl ikK 2 .t -K 4 ) I j I . (Ill) 



\mj 




If we restrict to j = jji(r, q, k), and m = *2j, then, acting on (r, q) G A4r,o, the 
only term of non-negative order is —K\, which can be expressed as K\r = d^G^r. 
Here 

G k r{x) = J g k (x - y)r(y)dA(y), (112) 



r 



is an operator of order —1. We use T^(k)(r, q) to denote this operator restricted 
to this subspace of data. 

The hybrid system of integral operators [46] is 



The range of ^dvT^k) is contained in the space of functions on T with mean 
zero on every component. Proposition Q] holds in the form version as well. 

Suppose now that k £ F + , and Q + (k)(r, q) = 0, with (r, q) G Mvfl \ {0}, 
and let (£ ■ , rj + ) be the solution to the THME(A;) defined by this data. The fact that 
Q + (k)(r,q) = 0, implies that 

d^+t = and V+ \ Tr = 0; (114) 

the second condition implies that dr£ +t = 0, as well. If the cohomology class 
[£+t]r = 0, then (£ + )t vanishes and Theorem [8]implies that (^ + , rj + ) is identically 
zero. When F is simply connected, _£fj R (r) = 0, and this proves Theorem |7J 
written in terms of forms. 
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Theorem 12. If Y is simply connected, then F + is disjoint from the closed upper 
half plane. Thus, the integral equation 




(115) 



provides a unique solution to the scattering problem from a perfect conductor for 
any k in the closed upper half plane. Here, 

f = G (d* v tf), ikhdA = d v if = ikrf* \ Tr , (116) 

where £J. n is the tangential component of an incoming electric field, and rj m |Yr > 
the normal component of the incoming magnetic field. 

When applying our method in the non-simply connected case, the following 
result is useful. 

Proposition 3. Suppose that k (£ E + U F + , and let tp 6 W 1 (r). The unique 
outgoing solution to the THME(k) with £ +t = ip is defined by data (r, q,j H ) with 

Proof. As k ^ E+, the proof of Theorem ITTI produces a solution, (^^y^ rj N+ ), 
to the THME(/c) with vanishing normal components and [£/v+t]r = [ip]r- The 
potentials corresponding to this solution take the form (ro,qo, j H ) with j H ^ 0. 
The condition [£jv +t ]r = [^]r shows that there is a function /, of mean zero on 
every component of T, satisfying 

Z N+t = Tp + drf. (H7) 

Since k ^ F + we can therefore solve the equation 

Q + (k)(r 1 ,q 1 ) = (0,God* r d r f). (118) 

With r} + ) the solution to the THME(/c) defined in SI by this data, we see that 

(£h+,Vh+) = (£ N +,Vn+) ~ (£+>»7+), (H9) 

satisfies 

dr£ H+t = d^H+t = 0, (120) 

and therefore £# +t G W 1 (r). This solution corresponds to the sources (rp — 
r i)9o — w i tri 3h ^ 0- Theorem [TOl then implies that £ H+ ^ 0. While 

it is not clear that [£# + t]r = [V ; ]r, it follows from (11201 ) and Theorem Q] that 
[Cn+tlr 7^ 0. Thus j H i — > %H+t i s an hijective linear mapping from W 1 (r) to 
itself, and therefore an isomorphism. □ 
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Remark 12. Let {ip 1 , . . . , ip 2g } be a basis for Ti 1 (T) . The proof of this proposition, 
along with that of Theorem [TT] show that, for k (£ E + U F + , we can effectively 
construct solutions { , 

V 1 +h) ■ ■ ■ » »7+h)}> in to the THME(fc), which 

satisfy 

i l +m = ^ l forl = l,...,2g. (121) 

Using Proposition [3j the system of equations, Q + (k), can again be used to 
solve the perfect conductor problem, at least for k ^ E + U F + . Let £™ be the 
tangential component of an incoming solution and set (/, h) as in d 1 1 6b - For any 
k ^ F+, with non-negative imaginary part, there is a unique solution to 

Q+(k)(r,q) = (f,h). (122) 

We let (£ + ,r? + ) be the solution to the THME(fc), defined in Q by this data. Let 
rj+) be the unique outgoing solution with 

£+t = C- (123) 

The tangential component of the difference tp = £ +t — £ +t belongs to W 1 (r). 
In the simply connected case this is zero, and therefore in this case we are done. 
In general, we can use Proposition [3] to find the unique solution , , f} + ), to the 
THME(fc) with £ +t = ip. By Theorem [T] the sum satisfies 

(€+,i7 + ) = + (€+,<)+), (124) 

and therefore solves the original boundary value problem. 

Remark 13. In our modification of T^~(k), we use —Go *2 as a "precondi- 
tioner," and to obtain a scalar equation. Other choices are possible, for example 
—Gi (*2^r + £) where / is another complex number and £ is an order zero operator 
mapping 1 -forms to functions. For numerical applications it may be important to 
find a good choice here. 

A question of considerable interest is to characterize the sets E + and F + . One 
might hope that E + is disjoint from the closed upper half plane. The set F + de- 
pends, to some extent on the choice of preconditioner. Indeed we can modify the 
preconditioner so that F + = E + . We let W be the L 2 -closure of {dru + d^(vdA) : 
u,v £ C°°(r)}, then, Theorem [3] shows that for k ^ E +1 W is a complement to 
the tangential ^-boundary values of the A;-Neumann fields: 

H k (n)t = {S +t : (£+,*?+) G W fc («)}. (125) 
To define the "optimal" preconditioner, we use the following lemma 
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Lemma 8. There is an analytic family of projection operators {Pk '■ k E C+ \ 
E + U {0}} satisfying 

• Pk \n k (Q) t = Id 

• Pk \w= 

Proof. The proof of Theorem [TT] gives an algorithm to construct a basis, 

{(&+(*0>»/W(*0) : l = l,...,2g} 

of TLk{^l): which depends analytically on k G C \ E + U {0}. Let Qo denote the 
orthogonal projection onto W 1 (r). Let {ip m : m = 1, . . . , 2g} be a fixed orthonor- 
mal basis for Tt l (T) . The matrix of the restriction Qo \n k (fi)t > w ^ tn res P ect to these 
bases, is given by 

A lm (k) = (e N+t (k)^ m )L H ry (126) 

This matrix is analytic and invertible in C + \ E + U {0}. The inverse transformation 
Rk : H l (T) — > TLk(Q)t is therefore also analytic. We define Pk to be 

P & (a) = R k Qo(a) (127) 

As Qo annihilates W, it is immediate that P& satisfies the conditions above. □ 

We can modify our hybrid system by letting 

As (ipPfc is a bounded, finite rank operator, the family Qf (k) is again Fredholm of 
second kind. Suppose that (k)(r, q) = for (r,q) ^ (0,0), and let (£ + ,r7 + ) 
be the solution of the THME(fc) defined by this data. We see that 

drH +t = and [d* r - d* r P k ]£ +t = 0. (129) 

If k ^ E + U {0}, then the boundary data of solutions in Hj.(Q) solve this system 
of equations. If there were another solution, then in fact we could find a 1-form 
a defined on T, which solves this equation, and P^a = 0. This means that a = 
dru + d^vdA, and 

d* r d r u = d r d* r vdA = 0, (130) 

which easily implies that a = 0. Thus for k G C + \ E + U {0}, TLk{^l)t is the 
complete set of solutions to the system of equations in d 1 29b and therefore the 
nullspace of (k) is trivial for k G C+ \ E + U {0}. Hence for this choice of 
preconditioner F + C E + . 

Additional care is required near k = 0, as the rank of the map £ + h-> [£ +t ]r 
drops at k = from 2g to g. 
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7.2 Low Frequency Behavior in the Non-simply Connected Case 

If T is not simply connected, then the space of solutions defined by data in Mr,o 
converges, as in the simply connected case, to the orthogonal complement of the 
span of the harmonic Dirichlet and Neumann fields. In this case we also need to 
consider what happens to the solutions defined by data from W 1 (r). Using this 
data we also obtain the harmonic Neumann fields. We consider the two types of 
data separately, beginning with that from Mr,o- 

Theorem [10] demonstrates that, for any k ^ E + with \mk > 0, we can solve 
the boundary value problem: 

d£ + = ikr 1+ d*r} + = -ik£ + 

= f i-n,(*3V+) = h, 

for arbitrary (/, h) in Air,o- Indeed as ^ E + , and the integral equations on V are 
of the second kind and analytic in /c, it follows that we can actually solve (11311 ) for 
k in an open neighborhood, V, of 0. We now discuss what happens to our solutions 
as k tends to 0, within a relatively compact subset of V. In particular, we would 
like to characterize exactly which harmonic fields arise as limits of fields of the 
form given in (l83T ). where (r, q) are obtained by solving (I33T ). j = jn(r, q, k) and 
m = * 2 j. 

To avoid confusion, we let (£ k , r) k ) denote the unique solution to (1131b . for a 
fixed (/, h) G Mr,o- As jji(r, q, k) is 0(k), it follows easily that, as k tends to 0, 
(,£kiVk) converges to 

£o = # Vo = d*$ m . (132) 

Theorem 13. The set of limits (£ , rj )for (/, h) G Mv,o is the orthogonal com- 
plement to the span of both the harmonic Dirichlet and Neumann fields. 

Proof. As the components decouple at k = it suffices to check each separately. 
Since the Hodge star-operator interchanges the solutions, as well as, the 1- and 
2-form Dirichlet/Neumann fields, we need only check the 1-form case. 

First we show that £ is orthogonal to the harmonic Dirichlet fields. These 
fields are of the form £ d = du, where u is a harmonic function, constant on each 
component of T. We observe that u = (|a;| _1 ) and £ = 0(|a;| -2 ), and this justi- 
fies the following integration by parts: 

(£o> £d)n = (£o> dr A u )r + (d*£ , u) n 
= (i n £o,u)r = 

The last equation follows as i n £ has mean zero over every component of T, and 
u is constant on each component. 
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We now turn to the Neumann fields. A Neumann field, £ n , satisfies £ n = 
0(\x\~ 2 ). In this instance we use the equation £ = d<p to conclude that, 

<£o> On = {dr A 0, £ n )r + {<P, d*Un 
= (<t>,in€ n )r = 

The last equality follows as i n £ n \r= 0, by definition. 

An outgoing harmonic 1-form is determined by it normal components along T, 
up to the addition of an arbitrary Neumann field. As the limit £ is required to have 
mean zero on every component of T, but is otherwise unrestricted, it follows that 
every outgoing harmonic 1-form, £, has a unique orthogonal decomposition as: 

€ = Co + €d + € n (135) 

One simply chooses £ d so that i n (£ — £ d ) has mean zero on every component of 
T. This then uniquely determines £ . Recalling that the Dirichlet and Neumann 
harmonic 1 -forms are themselves orthogonal, the Neumann component is then de- 
termined by orthogonally projecting (£ — £ d — £ ) ont o the Neumann fields. This 
completes the proof of the theorem. □ 

We now turn to data from W 1 (r). The first cohomology group of T splits into 
two disjoint subspaces, one is the image of the restriction map H dR (D) — > TL l (T) , 
the other the image of i7^ R (f2) —>■ Tt l (T). By the Mayer- Vietoris sequence, these 
restriction maps are injective, and so, by a small abuse of terminology, we may 
speak of 5^(0) and H^ R (D) as subspaces of iJ^ R (T), and write 

H 1 dR (T) = H 1 dR (n)®Hl K (D), (136) 

see ||29l . With this notation, H^ R (Q) is dual to the "A-cycles," shown in FigureQ] 
while H dR (D) is dual to the "B-cycles." 

The solutions to THME(O) are harmonic fields, which satisfy the decoupled 
equations: 

d£ ± = d*£± = and dr}± = d*r] ± = 0. (137) 

From these equations it is clear that £ ± and *3?7 ± are closed 1-forms and therefore 
define classes in their respective ff^-groups. We see that 

K +t ]r, [(*aV+)t]r G H^ R (n), and [£_ t ] r , [(*3V-)t]r G H dR (D). (138) 

The jump relations show that the solution of the THME(O) defined by the data 
(0, 0,j H ), satisfies 

[£+t - £-t] = -3h and K*3V+)t - i*3V-)t] = *2J H - (139) 
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Choose a basis of harmonic 1-forms, {ip : I = 1, . . . , g} for H^ R (Q) C 
fldn(r). Their Hod § e duals {^ 9+l = = i = 1, • • • ,9} 

are also harmonic 

and are a basis for H^D) C #dii( r )- We 

say that such a basis is adapted to 
the splitting of Hl R (r) in (fT36l We let r^) : Z = 1, . . . , 2g}, denote the 
solutions to the THME(O) defined by this data, with (r, q) = (0, 0). The relations 
in (11381 ) imply that image of each of the restriction maps £ + ^ [£+t]r> £ l— ► 
[£-t]r, spans a subspace of //j R (r) of dimension at most 5. On the other hand, 
the jump relations show that the differences, 



span all of H^ K (T). These relations, and analogous ones for the ^-components, 
along with d 139b . easily imply the following result: 

Proposition 4. The solutions (£±,77^), satisfy 

1. For I = 1, . . . , g, the restrictions [£+Jr span H^(VL) C H^(T), while the 
restrictions [(*3*7+)t]r = 0. 

2. For I = g + 1, . . . , 2g, the restrictions [(*3f7+)t]r span H^ R (Q) C H^(F), 
while the restrictions [£+Jr = 0. 

3. For I = g + 1, . . . , 2g, the restrictions [£'_Jr span H^ R (D) C H^ R (T), 
while the restrictions [(*3?/-)i]r = 0. 

4. For I = 1, ...,g, the restrictions [(*3T]L)t]r span H^ R (D) C H^ R (T), 
while the restrictions [£'_Jr = 0. 

We now recall the basis of /c-Neumann fields, JV* = {(^5v + (fc), r] l N+(^)) ' 
1 = 1,..., 2g}, constructed in the proof of Theorem [TU This is an analytic family 
in a neighborhood of 0, as it only requires the solvability of the normal equations. 
We now assume that we define these fields, using a basis, : 1 = 1,..., 2g}, 
of Tt 1 ^), which is adapted to the splitting in (1 1 36b - Proposition [4] shows that at 
k = the fields in Nk continue to span a 2<?-dimensional vector space of solutions 
to the THME(0), and that {£ l N+ (0) : I = 1, . . . , g} are a basis for the space of 
outgoing, harmonic 1-forms, with vanishing normal component along Ml. Note 
that for / = 1, . . . , g we have 



These fields are outgoing, harmonic 1-forms, with vanishing normal components 
along Ml, hence there must be constants {a%, . . . , a g } so that 

9 



(140) 



i„$?(0) = 0and (0)] r = 0. 



(141) 




(142) 



m=l 



40 



The second equation in (11411 ) implies that all the coefficients are zero. This proves 
that the basis Nf* reduces at k = 0, to a basis of the form: 



{(^ + (0),0), (0,7^(0)): l = l,...,g}. 



(143) 



We summarize these results in a theorem. 

Theorem 14. There is an open neighborhood U ofO £ C, and 2g analytic families 
of outgoing solutions {(£ l N + (k), ri l N+ (k)) : I = 1, . . . ,2g} to the THME(k), 
which, for each k £ U are a basis for the k-Neumann fields 7ik{£l)- At k = the 
and -q-components decouple and satisfy (1 143b - 

Theorems [13] and [14] give a clear picture of the behavior of the space of solu- 
tions to the THME(&;) in a neighborhood of zero, defined by the representation ([83l . 
They show that, in a reasonable sense, this representation does not suffer from low 
frequency breakdown. Let U C C be a neighborhood of zero, and {a(k) : k G U} 
be a continuous family of 1-forms defined on T. If a is orthogonal to W 1 (r) and 
dra(k) j k has a limit as k tends to zero, then it is clear that the hybrid system pro- 
vides a continuous family of solutions, in a neighborhood of zero, to the THME(A;) 
with £ +t (k) =<*(&). In a subsequent publication we will consider conditions on 
the projection of a(k) into Tt 1 ^), which are needed to conclude the existence of 
such a continuous family of solutions. 

8 The Normal Component Equations on the Unit Sphere 

In this and the following section we determine the exact form of the systems of 
Fredholm equations derived above for the special case of the unit sphere in R 3 . We 
make extensive usage of spherical harmonics and "vector" spherical harmonics, in 
the exterior form representation. As this is not standard, these formulae are derived 
in Appendix IA.5I The equations decouple, and very nicely illustrate the general 
properties described above. As the equations for the normal components are a bit 
simpler, we begin with them. 

The integral equations for the normal components of £ and r/, can be solved 
simply and explicitly when T is the unit sphere centered at 0. This reveals the close 
connection between our equations and the Mie-Debye solution. We are represent- 
ing £ and r\ in terms of the potentials a, a m , (j), and <E> m , with j a 1-form on S\ 
and m = *2j - The Debye sources r, q satisfy: 




(144) 
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If we assume that 



r = Y j ai m Yi m q = Y J h lmY l m 

lm lm 



(145) 



im 



then (11441) implies that 

- 1(1 + l)a im = ikai m and - 1(1 + l)(3i m = ikb im . 



(146) 



Suppose that the normal components of £ and rj are represented in terms of spher- 
ical harmonics by 



i£ = Yl ClmY i m and i n*z r n = Yl d ^ Y i m - 



(147) 



im 



lm 



Using the results of Propositions [7] and [U in the appendix, we see that the in- 
tegral equations in (I33T ) for the different spherical harmonic components decouple. 
The equation for the coefficient of the /m-component of i n £ becomes: 



r, V r 
Hm I l 



ik ( ai m G k 



(dslYD ■ 



dx 



(* 2 d s! Yn ■ dx\) 



+ PlmGk 

dG k (a lm Yn 
Msjir) • dx] - p lm G k [(d sl Yn ■ dx 
Using the results of these propositions we see that this equation reduces to: 



*3<^ ( OtlmGk 



(148) 



cim = ika 



im 



ikl(l + 1) 



21 + 1 



ik 2 a lm j l (k)d k h ( p(k) + ikl(l + l)a lrn jj(*0fy (1) (*0- (149) 
Standard recurrence relations for the spherical Bessel functions imply that 

k [j^h^k] - j^h^k)] 



21 + 1 



m ^ (k) -+ji(k)d k h\ 1 \k) + d k j l (k)h\ 1 \k), (150) 



k 
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see lfT3l . Using this identity and the relations in (11461) . we obtain 

cim = ai m khf\k) (iji(k) + ikj[{k) + kji{k)) . (151) 

We define the function 

m n (k,l) = khf\k) ((» + fc)iz(fc) +iA$(fc)) . (152) 

An essentially identical sequence of steps leads to the relations: 

dim = -m n (k,l)bi m - (153) 

The diagonal entries of the block diagonal matrix we need to invert to solve the nor- 
mal component problem for the unit sphere, at frequency k, are simply {m n (k, I) : 
1 = 1,2,...}. 

Remark 14. The connection to classical Debye theory is now easy to establish. If 
we expand the Debye potentials u, v in (0 as 



v{r,e,4>) = Ys^hl (kr)Y t 



l,m 



Lm 



then a straightforward calculation E21 shows that (in terms of the normal compo- 
nents) 

1 

aim ~ J(fTT) Clm 

Thus, our generalized Debye sources, defined only on the surface, are analogous 
(but not equivalent) to the restrictions of the Debye potentials to the sphere. The 
classical Debye approach requires that the potentials themselves (defined in M 3 \ T) 
be expanded in surface harmonics, preventing the approach from being extensible 
to arbitrary geometry. 

We apply the Wronskian identity, 

ji(k)dhf\k) - d k ji{k)hf\k) = p, (154) 
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to find that: 



m n (k,l) = (l + ^iK^/ip^^ + ^iK^^^W + ^iK^^W 



(!)/ 



,(!)/ 



(155) 

Applying standard asymptotic formulae for ^ and to this representation shows 
that, for a fixed k with non-negative real part, we have: 



(156) 



This agrees with the fact that the integral equation for i n £ is of the form ^ + K(k), 
where K is compact. It is also the case that 



m n (0,l) 



1 



2 2Z + 1' 



(157) 



which shows that these equations do not exhibit low frequency breakdown. For 
integral I, and k along the real axis we have 



m n (k,l) ~ 1 + 0{k~ 1 ), 
as Im k tends to infinity we have: 

m n {k,l) ~ l + CKAT 1 ). 



(158) 



(159) 



Figure [2]X) shows plots of {|m n (/c, Z)| : = 1, 10, 100}. The condition number 
increases with the frequency, k. 




20 40 60 80 100 120 140 160 180 200 



Figure 2: Plots of |m n (l,Z)|, |m n (10,/)|, |m n (100, Z)|. 
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For fixed I the solutions of m n (k, I) = have Irak < 0. If Z is fixed then the 
imaginary parts of the roots of m n (k, I) = 0, decrease in proportion to minus the 
log of the real part, 

1m k oc --logRefc. (160) 
The first 50 zeros of m n (k, 1), rn n (k, 5), and m n (k, 7) are shown in Figure[3] 



20 40 60 



100 120 140 



(a) Zeros of m n (k, 1). (b) Zeros of m n (k, 5). (c) Zeros of m n (k, 7). 

Figure 3: Graphs of the first 50 zeros of the multipliers m n (k, I) for I = 1,5, 7. 

Using the three term recurrence relations for spherical Bessel functions, {z{\, 
we easily obtain the functional equation: 



Zl (-k) = (-i)^(£0; 

using this identity and (11521 ) it is not difficult to show that 



m n (-k,l) = m n (k,l). 



The function can be factored as 



hf\k)= Pl {k) e 



lk 



k l+l 



(161) 



(162) 



(163) 



where pi is a polynomial of degree /. Thus the multiplier takes the form 

lk ({i + m^ + ikjiiky 



m n (k,l) = Pl (k)e l 



k l 



Pi(k)e l 



k ( { k - il)ji(k) +ikji_ 1 (k) 
k l 



(164) 



It is easy to see that the numerator has a zero of order I at k = 0, and therefore the 
quotient is regular and non- vanishing there. The polynomial contributes I roots; the 
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symmetry, (1 162b . shows that, when I is odd, one root lies on the negative imaginary 
axis. A plot showing the roots with smallest imaginary part, and positive real part 
is shown in Figure [H 



100 200 300 400 500 600 700 800 900 1.000 



Figure 4: The solutions of m n {k, I) = with smallest modulus, and positive real 
part, for I = 1 to 1000. The real part increases monotonely with I. 

Remark 15. Theorem [lOl shows that the frequencies k, with Imk > 0, for which 
equation (133T ) has a non-trivial null-space coincides with the eigenvalues of the 
interior boundary value problem for Maxwell's equations defined by 

£_ fr= i n V- It ■ (165) 

Of course there are no eigenvalues, or resonances with Im k > 0. Calculations like 
those above, though simpler, show that the boundary condition holds for the vector 
spherical harmonics of order I provided: 

(i + k)ji(k) + ikj[(k) = with k ^ 0. (166) 

From (1164 1 ), we see that the left hand side of (11661 ) is a factor of m n (k, I). Thus the 
eigenvalues of the interior problem are a subset of the resonances of the exterior 
problem. These eigenvalues are the "non-physical" interior resonances connected 
with our representation (I83T ) of £ and r] in terms of potentials. They are familiar 
from the EFIE and MFIE representation, but shifted to the lower half plane, where 
they do no serious harm. The roots of hf: (k) are known to be related to scattering 
resonances for scattering off of a conducting sphere. 



9 The Hybrid Equations on the Unit Sphere 

To find the precise form of the hybrid operator, Q + (k), on the unit sphere, we only 
need to work out —God* 2 *2 T^(k). The normal equation is given by (1153b - As in 
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the previous section, we represent £ and 77 in terms of the potentials a, a m , (ft, and 
$ m , with j = j R (r, q, k) a 1-form on 5f and m = *2j. The generalized Debye 
sources r, q satisfy (11451 ) and (11461) . The £-field is given, in terms of the potentials 

by 

£ = [ikG k j ■ dx - dG k r - * 3 dG k * 2 j ■ dx}. (167) 

Using the expressions for r and j in terms of spherical, resp. vector spherical 
harmonics, we see that the tangential components of £ , modulo ker d% 2 , are given 
by 



£ +t mod ker d* S 2 



l, in 



-k 2 ai r , 



2I + , j [(* + l)ji-i(k)hf\(k) + lj l+ ih\^(k)] 



-ikatmjtWhW \k) + ika^jiikWik) + kd k hY\k)} 
Using (11461 ) and the identity, 



,(1), 



we easily obtain that 



Go4 f T+(fc) 



S^~^ j 'I In: ' / 



(169) 



9/ <^\2Z + 1 



- ikl(l + l)ji(k)h^\k) + fe 2 Ji(fe)[/», (1) (fc) + W, (1, (*0] (170) 



+ 



2/ + 1 



As with the system of normal equations, the hybrid equations are decoupled, pro- 
viding one equation for the coefficients of r and one for the coefficients of q. The 
only term on the right hand side of d 170b that is not 0(Z _1 ) is 



-iklil + lMk^ik) 



-1 



+ o(r 1 ) J 



21 + 1 4 

in agreement with (|48T ). We use the identity 

d k {khf\k)) = khf\{k) - lhf\k) 



(171) 



(172) 
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to remove the derivative from (11701 ). and define the multiplier for the tangential 
equation: 



m t (k,l) 



-k 



il(l + l)j,(fc)fy (1) (A;) - kji(k) [kh^^k) - Ih^'ik)} 



(i) 



21 + 1 

ik 2 [(l + l)j/-i(fc)/ig 1 (fc) + Ij^Wh^jk)] 
21 + 1 



(173) 



This multiplier behaves much like the multiplier m n (k, I) found for the normal 
equations. For fixed I its roots, as a function of k, lie in the lower half plane. The 
plots in Figure [5] show contours of log \mt(k, l)\ for I = 1, 10, 20. The x-axis is 
shown as a black horizontal line. They clearly show that the zeros lie in the lower 
half plane, and indicate the moderate behavior of the multiplier in the upper half 
plane. The plots in Figure [6] show the \m t (k, l)\ for I between 1 and 20, and 20 
and 200, respectively, for and k = 1, 10 and 100. It should be recalled that there is 
a certain amount of arbitrariness in the definition of this multiplier, resulting from 
the arbitrariness in the choice of preconditioner, God% 2 in the present instance. 

If the incoming tangential data are given by 



Lm 



Pi i 



* 2 d sl Y{> 



(here we used the normalized basis elements) then 

*2d S 2^ * 2 d S 2 * 2 d S 2Y t 



(174) 



*3 ^7 



ik 



Lm 



Hr, 



iky/l(l + 1) 



l.m 



(175) 



To find the tangential data for the hybrid system we apply God% 2 to £i n , ob- 



taming: 



God^t 



V #Tl) a y,. 



Lm 



2I + 1 



(176) 



For the unit sphere, the hybrid equations are therefore: 

m n (k,l)k m 
mt(k,l)a lm 



ik 



(177) 



21 + 1 



Ph. 
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These equations display a mild sort of low frequency breakdown, in that the co- 
efficients of the normal data, + l)<7zm}> must be uniformly 0(lo) in order 
for this system of equations to be stable. Of course the incoming data 77™) is 
assumed to be a solution of the THME(A:), so these estimates should automatically 
hold. Indeed if T7" 1 is given, then there is no need to differentiate and divide by 
k to find the data for the normal equation. 




5 ID 15 20 25 30 5 10 15 20 25 



(a) Contour plot of of log \m t (k, 1) | . (b) Contour plot of of log \m t (k, 10) | . 




5 10 15 20 25 30 



(c) Contour plot of of log \mt(k, 25) |. 

Figure 5: Plots of log \rnt(k, 1)|, log \mt(k, 10) |, log \m t (k, 25)|. The black hori- 
zontal line indicates the x-axis. The zeros are located near the deep blue dots. 

The use of Go in the preconditioner also leads to growth in the multiplier 
mt(k,l) are k increases for fixed Figure [(Jc) show \mt(k, 1)|, \mt(k, 10)|, and 
\mt(k, 20)|, for real k 6 [0, 100]. In the interval < k < I these functions oscil- 
late around a small non-zero value. When k exceeds / these functions show linear 
growth. Replacing Go, with something like Gji fc i should fix this problem. 
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(a) Small / 



(b) Large I 




(c) Fixed I 

Figure 6: Plots of \m t (k, 1)|, for a large range of I and k = 1, 10, 100, and fixed 
/ = 1,10,20 with k G [0,100]. 

10 Conclusions 

In this paper, we have developed a new representation for solutions of the time 
harmonic Maxwell equations, exterior to closed surfaces, based on two scalar den- 
sities. In the zero frequency limit, these densities are uncoupled and correspond 
to electric and magnetic charge. At non-zero frequency, however, they do not cor- 
respond directly to physical variables. They are simply used to construct electric 
and magnetic currents, after which the classical scalar and vector potentials and 
anti-potentials are employed (in the usual Lorenz gauge). Because of the close 
connection to the Lorenz-Debye-Mie formalism when the analysis is restricted to 
the unit sphere, we refer to our unknowns as generalized Debye sources. The natu- 
ral boundary data for our unknowns are the normal components of the electric and 
magnetic field and we have provided a detailed uniqueness theory for this bound- 
ary value problem for boundary surfaces of arbitrary genus (Theorems |U [5j |9]). 
In the course of this analysis, we have given a new proof of the existence (in the 
non-simply connected case) of families of nontrivial solutions with zero boundary 
data, which we refer to as fc-Neumann fields. They generalize, to non-zero wave 
numbers, the classical harmonic Neumann fields (Theorem ITTb. 
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We have also introduced a new Fredholm integral equation of the second kind 
for scattering from a perfect electrical conductor and have shown that it is invert- 
ible (in the simply connected case) for all wave numbers k in the closed upper 
half plane. There is a natural extension of the approach to the case of a dielectric 
interface, which will be reported at a later date. 

The work begun here gives rise to a new set of analytic and computational 
issues. In order to use the Debye sources as unknowns, one needs an efficient and 
accurate method for inverting the surface Laplacian. For surfaces T of genus g > 0, 
we also need to be able to efficiently construct a basis for the harmonic forms 
H l (T) . Finally, additional work is required to extend our approach to open surfaces 
(see, for example, [14]). These arise as common and important idealizations in the 
analysis of thin plates, cylindrical conductors and metallized surface patches in 
radar, medical imaging, chip design and remote sensing applications. 



Appendix 

A Exterior Forms, Maxwell Equations and Vector Spher- 
ical Harmonics 

In the traditional approach to electricity and magnetism Maxwell's equations are 
expressed in terms of relationships between four vector fields E, D and B,H 
defined on M 3 x R : 

dD dB 
— = cV x H - AttJ — = -cV x E 
dt dt (178) 

V • D = 4vrp V • B = 0; 

c is the speed of light. Here J is the current density and p is the charge density, 
they satisfy the conservation of charge: 

| = -V.J. (179) 

The differential symmetries of this system of equations are rooted in the exactness 
of the sequence: 

C°°{U) C°°(U] TR 3 ) ^ C°°(U; TR 3 ) ^ C°°{U), (180) 

here U C M 3 is an open set, and C°°(U; TM 3 ) are the smooth vector fields defined 
in U : 

C°°([/;TR 3 ) = {ai(x)d xi +a 2 (x)d X2 +a 3 (x)d X3 : ai ,a 2 ,a 3 £ C°°(U)}. (181) 
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The exactness of the sequence is equivalent to the classical identities V x V = 0, 
and V • Vx = 0. 

While this representation is traditional, physically and geometrically it makes 
more sense to regard Maxwell's equations as a relationship amongst differential, or 
exterior forms. In the second part of this paper we usually work with the fields E 
and H. It turns out to be convenient to use a 1-form to represent E and a 2-form 
to represent H. We use the correspondences 

H = h\d Xl + fi2d X2 + h 3 d X3 <-> h\dx 2 A dx 3 + h 2 dx 3 A dx\ + h 3 dx\ A dx2 = rj 
E = e\d Xl + e 2 d X2 + e 3 d X3 <-» e\dx\ + e 2 dx 2 + e 3 dx 3 = £ 

(182) 

Under this correspondence £ is the metric dual of E and *rj (* is the Hodge star- 
operator defined by the metric on M 3 ) is the metric dual of H. 

It is natural to think of the electric field as a 1-form, for the electric potential 
difference is then obtained by integrating this 1-form: 

<\>p ~ <t>Q = Jt (183) 

7 

with 7 a path from P to Q. Similarly, it is reasonable to think of the magnetic field 
as a 2-form, for the flux of H through a surface £ is then obtained by integrating: 



Flux of H through S = 




(184) 



s 



While these are the most basic measurements associated to electric and magnetic 
fields, there are times when it is natural to integrate the .E-field over a surface, or 
the .tf -field over a curve. This is done, in the form language, by using the Hodge- 
star and interior product operations, i v , introduced below. A detailed exposition 
of this approach to Maxwell's equations can be found in O. 

For the sake of completeness we recall the definition of d on forms defined on 
M 3 : 

0- forms A°IR 3 : da = d Xl adx\ + d X2 adx 2 + d X3 adxs 

1- forms A 1 ]^ 3 : d{a\dx\ + a 2 dx 2 + 03^3) = (d xi a 2 — d Xl ai)dx± A dx 2 
+ (d X3 ai - d Xl a 3 )dx3 A dx\ + (d^as - d X:j a 2 )dx 2 A dx 3 

2- forms A 2 M 3 : d{a\dx\ A dx 2 + a 2 dx 3 A dx\ + a 3 dx 2 A dx 3 ) = 
{d Xl ai + d X2 a 2 + ^303)^x1 A dx 2 A dx 3 

3- forms A 3 M 3 : d{adx\ A dx 2 A dx 3 ) = 0. 

(185) 



52 



The sequence (11801 ) becomes: 

C°°(U) C°°(U; AV) C°°(U; A 2 M 3 ) C°°(U; A 3 M 3 ). (186) 
All of the classical differential relations are simply d 2 = 0. 

A.l Exterior forms on a manifold 

Generally we can define the smooth exterior p-forms on an n-dimensional mani- 
fold M, C°°(M; A P T*M), as sections of the vector bundle A P T*M. The exterior 
derivative is a canonical map 

d : C°°(M; A P T*M) -» C°°(M; A P+1 T*M). (187) 

If (xi, . . . , x n ) are local coordinates, then a p-form can be expressed as 

a = ai(x)dxi 1 A • • • A dx ip , aj G C°°. (188) 

Here T p is the set of increasing p-multi-indices l<ii<---<z p <n. In these 
local coordinates, the exterior derivative of a function f(x) is defined to be 

n 

df = ^d X] f(x)dxj, (189) 
i=i 

and of a p-form 

da = ^2 dai(x) A dx^ A • • • A (190) 

The remarkable fact is that this is invariantly defined; though it is really nothing 
more than the chain rule. The fact that mixed partial derivatives commute easily 
applies to show that d 2 = 0. If a and /3 are forms, then we have the Leibniz 
Formula: 

d(a A /3) = (da) A (3 + (-l) deg "a A d(3. (191) 

If a is a /c-form defined in an open subset, U of M then we can integrate it 
over any smooth, oriented compact submanifold S fc CC U, of dimension k, with 
or without boundary. We denote this pairing by 




(192) 



s 
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If a is an exact form, that is, at = d(3, and £ is a smooth, oriented submanifold 
with or without boundary, then Stokes' theorem states that 

J d(3 = J /3. (193) 

s as 

The boundary must be given the induced orientation. Note, in particular, that if 
<9£ = 0, then the integral of d(3 over £ vanishes. 

There is a second natural operation on exterior forms that satisfies a Leibniz 
formula. If v is a vector field, then the interior product of v with a fc-form, u, is a 
(k — l)-form, i v u, defined by: 

i v u{v 2 , ...,v k ) = u(v, v 2 ,..., v k ). (194) 

If lo and T) are exterior forms, then 

i v [co A rj\ = [i v cj] Ar]+ (-l) dc ^io A [i v rj]. (195) 

Using forms simplifies calculations considerably because forms can be auto- 
matically integrated over submanifolds of the "correct" dimension, keep track of 
orientation, and all the differential relationships follow from the fact that d 2 = 0. 
Moreover, Stokes' theorem subsumes all the classical integration by parts formulas 
is one simple package. 



A.2 Hodge star-operator 

To write the Maxwell equations we need one further operation, called the Hodge 
star-operator. This operation can be defined on an oriented Riemannian manifold. 
Suppose that ui, . . . ,to n is an local orthonormal basis of one forms, and 

dV = ui A ■ ■ ■ Au n , (196) 

defines the orientation. Ifl<ii<---<i p <n, and j\ < ■ ■ ■ < j n - p are 
complementary indices, then we define 

★ [u h A • • • A uj ip ] = {-lyuj, A • • • A uj jn _ p , (197) 

with e = or 1 , chosen so that 

u h A ■ ■ ■ A u ip A A • • • A u ip ] = dV. (198) 
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For example, the Hodge operator is defined on the standard orthonormal basis 
of exterior forms for M 3 by setting: 

★1 = dx\ A dx2 A dx3 
■kdxi = dxo A dx% * dxo = dx* A dx\ * dx* = dx\ A dxo 

(199) 

■kdx\ A dx2 = dx% * dx% A dx\ = dx% * dx-i A dx 3 = dx\ 
*dx\ A dx2 A dxa = 1 

From these formula it is clear that, in 3-dimensions, -k 2 = I. Notice that applying 
the ^-operator exchanges the two correspondences between vector fields and forms 
in (TT821) . 

The star-operator is simply related to the metric: If a, (3 are real p-forms, then 
[aA*/3] x = (a,f3) x dV x , (200) 
where (•, •) is the (real) inner product defined by the metric on p-forms. Generally, 

* 2 = (_1)P("-P) s (201) 

on p-forms defined on an n-dimensional manifold. Using this observation and (12001 ) 
we easily show that * is a pointwise isometry: 

(a,0) x = (*a,*/3) x (202) 

A fundamental role of the Hodge ^-operator is to define a Hilbert space inner 
product on forms. If a and (3 are (possibly complex) forms of the same degree 
defined in U, then a A */3 is a n-form, which can therefore be integrated: 



(a, (3) = J a A*/3. 



U (203) 
= J (ot,(3) x dV(x). 

v 

We assume that (•, •) is extended to define an Hermitian inner product on complex 
valued forms. The extended metric continues to satisfy (1202b . 

A.3 Adjoints, Integration-by-parts and the Hodge Theorem 

On an n-dimensional manifold the expression for the formal adjoint, with respect 
to the pairing in (12031 ). of the <i-operator, acting on ap-form /3, is: 

I — ★ d * 3 if n is even 

d*3 = < ^ (204) 

|(-l) p *d*/3 if n is odd. 
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Let G be a bounded domain with a smooth boundary. Let r be a function that 
is negative in G and vanishes on bG. Suppose moreover that {dr, dr) x = 1 for 
x G bG, and n is the outward pointing unit normal along bG. The basic integration 
by parts formulae for d and d* can be expressed in terms of this inner product by 



(da, (3) G = J {dr A a, (3) x dS(x) + (a, d*(3) G 

bG (205) 

(d*a,/3) G = - {i n a,(3) x dS{x) + (a,df3) G . 



bG 

Here we use dS to denote surface measure on bG. It is important to recall that, with 
respect to the pointwise inner product, 

{dr Aa,(3) x = (a,i n [3) x . (206) 

The (positive) Laplace operator, acting on any form degree is given by formula 

dd* + d*d. (207) 

In i? 3 this would give —{d Xl + d X2 + d X3 ). To avoid confusion with standard usage 
in E&M, we use A to denote the negative operator — {dd* +d*d). If M is a compact 
manifold without boundary, then the de Rham cohomology groups are defined, for 

< k < dim M, as 

= kev{d : C°°{M; A k T*M) -» C°°{M; A k+1 T*M)}/dC°°{M; A fc_1 T*M). 

(208) 

It is a classical theorem that these abelian group are topological invariants, see [29]. 
We let H k {M) denote the nullspace of the Laplacian acting on C°°{M; A k T*M). 
Stokes' theorem shows that 

(Aw,w) = (du,duj) + (d*Lu,d*uj). (209) 

Thus TL k {M) ~ ker d n ker d*. The Hodge theorem states that 

Theorem 15 (Hodge). If{M, g) is a compact Riemannian manifold, without bound- 
ary, then, for each < k < dim M, 

H k {M) ~ H^ R {M). (210) 

and, as an I? -orthogonal direct sum, we have: 

C°°{M;A k T*M) = dC°°{M; A fc_1 T*M) © d*C°°{M; A k+1 T*M) © 7i k {M). 

(211) 
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As A is an elliptic operator, the Hodge theorem has the following very useful 
corollary: 

Corollary 6. If(M, g) is a compact Riemannian manifold, without boundary, then 
for each < k < dimM, the dim H^ R (M) is finite. 

A.4 Maxwell's Equations in terms of exterior forms 

With these preliminaries we can state the correspondences between the differential 
operators, V, Vx, and V- and the corresponding objects acting on forms. For a 
scalar function <f>, V</> corresponds to dtp. An elementary calculation shows that if 
E <-> £, a 1-form, then V x E <-> d£, and V • E <-> d*£. Moreover with H <-> -q, 
a 2-form, we have V x H <-> d*rj, and V • H <-> drj. The operator d* also acts on 
3-forms. 

If we let E <-> £, ii" <-> 77, and J <-> j (a. 1-form) as in (1 1 82b . then Maxwell's 
equations in a vacuum become: 

9 t A* A ^ At 

— = cd r\ - 4ttj — = -cd£ 

at at (212) 

d*Z = 4Trp d V = -£ = -d*j. 

If £ and r] are time harmonic with time dependence e~ ltuJ , then in the absence of 
sources, we easily derive the Helmholtz equations: 

c 2 A£ + u; 2 £ = c 2 Arj + uj 2 rj = 0. (213) 

We let Dei 3 denote a bounded set with smooth boundary and let 

Q = R 3 \D r = bD. (214) 

In this paper D is usually taken to be a perfect conductor, lying in a bounded 
domain with smooth boundary, and a dielectric. We assume that e is the electrical 
permittivity, // is the magnetic permeabilty and a the electrical conductivity of O. 
As above, we identify the Enfield with a 1-form, S : 

e\d Xl + e 2 d X2 + e 3 d X3 = E <-> 3 = eidxi + e 2 cfa;2 + e 3 dxs, (215) 

and / / with a 2-form, AT : 

hid Xl +h,2d X2 + 1138x3 = H <-> AT = h\dx2 f\dx3 + h 2 dx3 Acfai + ^3^1 Adx 2 - 

(216) 
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In terms of exterior forms, we set 



l(x,t) 



uJ 



£(x)e 



-iuit 



uje + la 

the time harmonic Maxwell equations become: 

d£ = ikq d*£ = 
d*rj = —ik£ dr\ = 0. 



N(x, t) = fj, 2 rj(x)e 



-iwt . 



(217) 



(218) 



Here k is the square root of fi(eio 2 + iau>), with non-negative imaginary part. 

With this choice of correspondence between the vector and form representa- 
tions, we can write the Maxwell equations in a very succinct and symmetric form: 



(d + d*)(£ + ri) =ikA(£ + r 1 y, 
here A is the operation defined on forms by 

A(a) = {-l) dcga a. 
Simple calculations shows that 

(d + d*)A = -A(d + <f) and A 2 = /, 



(219) 

(220) 

(221) 
(222) 



implying that 

(d + d* — ikA) 2 = -(A + k 2 ). 

Thus, acting on forms, the operator A + k 2 has a local square root. Or, put differ- 
ently, d + d* — ikA is an operator of Dirac-type, see JH. Indeed, we could write 
the vacuum Maxwell equations in the form 



noting that 



[c{d + d*) + Ad t ](Z + n) = 4ircp, 
[c(d + d*) + Ad t } 2 = d?-c 2 A. 



(223) 
(224) 



A.5 Vector Spherical Harmonics 

The domains defined in M 3 as complements of a round sphere are very important 
in applications. They also provide a context where the integral equations defined in 
the earlier sections can be diagonalized and solved explicitly in terms of classical 
special functions. In this appendix we give a brief treatment of the theory of vector 
spherical harmonics in the exterior form representation. A classical treatment is 
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given in sections 9.6-7 of lfT3l . But for the classical theory of (scalar) spherical 
harmonics (which can also be found in Jackson), our discussion is essentially self 
contained. 

We begin with the relationship between the (negative) Laplace operators in M 3 
and on the unit sphere Sf CK 3 . Recall that on any Riemannian manifold, (X, g), 
the Laplace operator on /c-forms is given by A£ = —{d*d + dd*), where * is 
defined by g. The following is simply the usual change of variables formula for 
spherical polar coordinates. 

Proposition 5. Let r 2 = x\ + x\ + x\, the scalar Laplace operator on M. 3 can be 
expressed as 

• is 1 p 20 , 1 A Sf , a Xid Xl + x 2 d X2 + x 3 d X3 

An = -r:O r r o r + -^A„ , where o r = - -. (225) 

We also need to relate the Laplacians on 1 -forms. 
Proposition 6. Let a be a 1-form on M 3 such that id r oc = 0, then 

Afa = -^Apa + L r a + (d* S 2<xj dr, (226) 

where 

2 

L r cx. = i r 2Q r d(i r -2Q r da) H — ^o.. (227) 

This formula follows by a calculation using a local co-frame field. If we ex- 
press a = J2 ctjdxj, then 

3 

L r ol = ^ ( d r a j ) dxj . (228) 

While ( 12261 ) is a good deal more complicated than ( 12251 ), it allows for a careful 

s 2 

analysis of the eigenforms of A x x . 

s 2 

We begin with the standard description of the eigenspaces of A 1 . Let £? de- 
note the linear space of scalar eigenf unctions on Sf satisfying: 

Af / = -1(1 + 1)/. (229) 
These spaces are represented in terms of classical spherical harmonics by 

£f = span{^ m : m = -/,..., 1} (230) 
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s 2 

A basis of eigenforms of A 2 L , with eigenvalue 1(1 + 1), is given by 

{* 2 Y™ = Y l m dA: m = -l...l}. 

The fact that i^ 1 (S' 2 ;lR) = and the Hodge theorem imply that the 1 -forms on Sf 
are the L 2 -orthogonal direct sum 

C^iSfiA 1 ) = d S 2C co (Sl;A )ed* s2 C co (Sj;A 2 ). (231) 

This observation coupled with the fact that d and d* commute with A imply that 

s 2 

the eigenspaces of A-^ 1 are given, for Z € N, by 

£} = span [{d S 2 : m = -I, . . . ,1} ffi fadcqYi 71 : m = -/,..., /}] , (232) 

with eigenvalue — /(/ + 1). From this representation, the classical orthogonality 
relations are quite easy: 

(d sf Yr, d sl Yf) = (d* s2 d s! Yr, Yf) = 5 u >5 mm 'l(l + 1) 

(* 2 d s! Yr, *2d sl Yp') = (dstYT, d s! Y t r') = 5 w 5 mm 'l{l + 1) (233) 

(d s2 Yr,*2d s! Yf) = (Yr^^Yf) = 0. 

The second line from the fact that *2 is an orthogonal transformation, and the last 
relation follows from Stokes theorem. 

The eigenforms ^dgiY™ are divergence free. If we extend them to M 3 so they 
annihilate d r , and express them in the form -kzdY™ = a\dx\ + a 2 dx 2 + a 3 dx 3 , 
where the coefficients are extended to be homogeneous of degree 0, then it follows 
from (12261 ) and the equation 

A? 3 (a 1 dx 1 +a 2 dx 2 +a 3 dx 3 ) = (A|f ai)dxi+(A|f a 2 )dx 2 +(Af a 3 )dx 3 (234) 

that is, where r = 1, we have 

A^oj = -1(1 + l)aj. (235) 

In other words the coefficients of ^dgzY™ lie in £®. These eigenforms correspond 
to the classical eigenfields of the form {r x VY™}. 

The members of the other family, {dg^Y™}, which corresponds to {r x (r X 
VY™)}, are not divergence free and their coefficients with respect to dxj lie in 
£®_ l © £i + i- These coefficients are easily found; if we think of Y™ as a homoge- 
neous function of degree on M 3 , then iQ r d^,Y[ n = 0, 

3 f)Y m 

d s 2 Yi m = £ ftT^i • (236) 
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In order to determine the action of the Green's function on the coefficients of these 
forms, we need to represent them in terms of spherical harmonics. Let 17™ = 
r l Y™ . This is a homogeneous harmonic polynomial of degree I. We see that 

FtTjm f)Y m t ■ 

^ = r l ^ + l^Ur. (237) 
oxj oxj r A 

We apply the Laplace operator to Xj U™ to obtain 

f)TT m 

Af( Xj un = 2^-, (238) 

and therefore 

x 3 UP = u% +r 2 u™J. (239) 

Here u^P are homogeneous harmonic polynomials of degrees I + 1 and / — 1, 
respectively. Once again applying the Laplace operator to this relation, we see that 

1 f)TT m 



and therefore 



J2 <?dx 3 = Urrdr - ^dlJ? . (241) 
j'=i 



Using the homogeneity we also see that 

~ 2Z + 1 



X>< = (242) 



Restricting to r = 1, gives 



W '-'= ^TiSr r - • (243) 

The functions on the right hand side belong £f_ Y and £^ +1 respectively. Employing 
these relations, we can work out the action of outgoing Green's function on £j. 
The outgoing Green's function for frequency k, with Im k > 0, is given by 

e ik\x-y\ 

g k (x,y) = — -. (244) 

4ir\x — y\ 

If r = \x\ > 1 and \y\ = 1, then we can expand g as 

5fc 



(x,y)=ikJ2jl(k)h ( i 1) (kr)P l , (245) 



z=o 
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here Pi is the orthogonal projection onto £®, and 

3i(z) = ]f[ z Jl + ^) and h?\z) = (z), 
see IPT3TI . When k = formula d245l ) reduces to 



(246) 



We let 



»(*.w) = E(2/ + l) r i+r 
G k f(x)= g k (x,y)f(y)dA(y). 



(247) 



(248) 



If q; is a 1-form on S^, then it has a unique extension to TM 3 f 5 2 that an- 
nihilates d r , which we denote by a ■ dx. The extended form has a well defined 
representation along S\ as 



a • dx = ajdxj. 

3=1 



(249) 



If we extend the coefficients to be homogeneous functions of degree zero, then 
iQ r a • dx = implies that 



3=1 



-a ■ dx. 



(250) 



which will prove useful below. 



Proposition 7. If \y\ = 1 and r = \x\ > 1, then, applied component-wise, the 
action of G k is given by: 



G k Yr = ikj l (k)h ( l l) (kr)Yf 



G k 
G k 



(* 2 d s? Yn ■ dx 

{d sl Yn-dx 



ikji(k)h\ 1] (fcr)(* 2 d 5 2^ w ) • dx 



ik 



Urdr 



dUp 
21 + 1 



[(I + (fc)/i[i\(fcr) + Ij^m^ikr)] 



(i) 



4-1 



^+l(fc)^ ( i ) 1 (fcr) 



(251) 
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On the right hand sides of (125 lb . Y™ is homogeneous of degree zero, as are the 
coefficients ofd^Y" 1 ■ dx, and ^d^Yj" 1 ■ dx. As above, UJ 71 is the homogeneous 
harmonic polynomial of degree /, defined by Yf"". 

Along the unit sphere the normal components are given by 

i 9r dG k Yr = i^mdkh^i^Yr 



ids**? 



dx 



ik 



1(1 + 1) 
21 + 1 



jl-i(k)hf\{k) - j l+1 (k)h\^(k) 



(i) 



Yi m - 
(252) 



Along the unit sphere the tangential components are given by 

[dG k Yr\ \ TS! = ikmhf\k)d sl Yr 
~G k K* 2 d s2 Yr) ■ dx]] \ TS 2= ikmhfXk) * 2 d sl Yr 



G k [(d sl YD • dx] 
ik 



(253) 



Its? 1 



We also need to compute the effect of *%d on these eigenforms. 
Proposition 8. Along the unit sphere we have: 



(i) 



d s ?Yr- 



id r *3 dG k 



{^ds^Y^) • dx = -ikl{l + l)j l (k)h ( l\k)Y l 



idr *3 dG k (dgaYD ■ dx 



0. 



(254) 



Finally we need to calculate the tangential components of these forms; the identity 
satisfied by a 2-form 

*3 « \ts*= *2[idr<x \ts*\i ( 255 ) 
along with (12501) facilitate these computations. 
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Proposition 9. Along the unit sphere we have: 

^ 3 dG k [(* 2 d s2 Yn ■ dx}] \ TS 2= -ikd S 2Yrji{k)[hf\k) + kd k hf\k)\ 



* 3 dG k [(d s! Yn ■ dx) 



ik * 2 d s 2Y{' 



i{i+2)j l + 1 {k)hfhk) 



TSJ ~ V 21 + 1 

_(/_!)(/ + i^^h^k) + m + m-^dkh^k) 

+ 1^+1^)8^^)] 

(256) 

Given (12351) . (I24TT) . (l243i (12421) . (|245T >. and (|253T > the formulse in these propo- 
sitions are elementary calculations, which follow from the fact that 

P t \ e o= l £ o and Pi \ £ o= Oifkj^l. (257) 
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